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Preface 


The  engineering  design  process  involves  a  mixture  of 
theortical  facts  and  experience.  This  experience  allows  many 
design  parameters  to  be  set  at  acceptable  levels  without 
detailed  exploration.  This  experience  factor  is  necessary 
since  any  design  problem  involves  almost  unlimited  possible 
configurations.  This  investigation  focused  on  developing 
anaylitical,  statistical  design  tools  to  help  specify  design 
parameters  for  a  robust  decoupled  controller,  an  area  with 
alimited  pool  of  engineering  experience.  While  significant 
progress  was  made  in  applying  these  statistical  tools  some 
judgements  were  still  necessary  on  such  items  as  ranges  of 
parameters,  and  scope  of  investigaiton.  I  myself  don't  have 
a  large  pool  of  experience  to  draw  upon.  This  factor  should 
be  kept  in  mind.  While  this  inexperience  may  have  reduced  my 
effectiveness  in  judgemental  decisions,  my  fresh  look  should 
be  of  benefit  especially  in  developing  unbasised  techniques. 

I  would  like  thank  my  family,  friends  and  church  for 
their  support  and  encouragement.  I  would  like  to  thank  my 
advisor  Major  David  G.  Robinson,  P2  and  my  committe  members 
Robert  A.  Calico,  Jr,  PhD,  and  Major  Kenneth  W.  Bauer,  PhD 
for  imbuing  me  with  knowledge  and  assisting  me  in  this 
effort . 


Jeb  Ewell  Brewer 
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Abstract 


This  investigation  developed  and  tested  statistical 
design  methods  for  configuring  a  robust  decoupled  controller 
on  a  lightly  damped  structure.  The  test  model  used  was  a 
lumped  mass  finite  element  representation  of  the  Charles 
Stark  Draper  Laboratory  Model  I  (CSDL-I).  The  decoupled 
control  methods  consisted  of  a  system  of  three  individual 
subcontrollers  designed  on  the  basis  of  a  subset  of  plant 
dynamics  with  control  authority  enacted  through  the  total 
system  of  structural  actuators.  Transforms  were  applied  to 
the  subcontrollers  to  insure  that  dynamic  coupling,  called 
observation  spillover  and  control  spillover,  didn't 
destabilize  the  global  system.  The  decoupled  controller 
system  was  designed  based  on  only  the  first  eight  of  the 
twelve  natural  frequencies  (modes)  for  the  CSDL-I.  The 
remaining  four  high  frequency  modes  were  modeled  as 
residuals . 

This  investigation  used  signal -to-noise  ratios, 
orthogonal  numerical  experiments,  and  linear  regression  to 
efficiently  probe  the  design  space  and  to  produce  a  robust 
control  system.  The  measure  of  system  performance  was  the 
structural  alignment  of  a  particular  node  referred  to  as 
1 ine-of-sight  (LOS).  This  LOS  requirement  simulated  space 
vehicle  requirements  associated  with  aligning  antennas,  LASER 
beams,  solar  arrays,  etc.  The  first  major  design  parameters 
were  the  elements  in  diagonal  state  weighting  matrices. 


These  state  weighting  matrices  were  used  in  Linear  Quadratic 
Regulator  and  Linear  Quadratic  Estimator  Theory  to  compute 
control  gains  for  the  subcontrollers.  The  second  major 
design  parameters  were  the  groupings  of  modes  within  the 
individual  subcontrollers.  A  heuristic  was  developed  to 
select  these  groupings  so  that  the  total  system  was  robust  to 
certain  noises.  The  noise  variables  investigated  included 
error  readings  in  position  sensors,  and  the  error  associated 
with  only  designing  the  system  based  on  eight  of  the  possible 
twelve  modes. 

It  was  concluded  that  the  techniques  were  efficient  in 
configuring  a  decoupled  controller  to  be  robust  to  sensor 
noise  and  modal  truncation  while  still  maintaining  adequate 
LOS  performance.  The  techniques  reduced  the  CPU  time 
required  from  a  possible  76,354  years  for  a  global  search  to 
only  16.5  hours.  It  was  also  found  that  the  robust 
controller  used  a  significant  amount  of  actuator  authority 
and  it  was  recommended  that  further  work  be  done  to  expand 
the  statistical  design  techniques  to  include  actuator 
authority  as  an  output  metric. 


STATISTICAL  TECHNIQUES  FOR  DESIGNING  A  DECOUPLED 
CONTROLLER  TO  BE  ROBUST  TO  MODEL  AND  SENSOR  NOISE 


I.  Introduction 


In  the  future  exploration  of  space,  larger  more  complex 
structures  will  be  placed  in  orbit.  Whether  these  structures 
will  be  manned  space  stations  or  unmanned  satellites  they 
most  probably  will  be  extremely  light-weight,  flexible 
strictures  requiring  accurate  structural  alignment.  High 
precision  structural  alignment  requirements  for  LASER  beams, 
solar  arrays,  antennas,  etc  will  necessitate  active 
structural  control  of  the  orbiting  vehicles.  Light-weight 
flexible  structures  possess  many  closely  spaced  low  frequency 
vibrational  modes  that  make  active  feedback  control 
computationally  burdensome.  In  an  effort  to  actively  control 
such  structures  while  minimizing  on-line  computational  loads 
two  design  approaches  have  been  used;  simplifying  the  dynamic 
model  representation,  and  using  decentralized  control 
methods.  Mercadal  (8:1-20)  provides  a  good  summary  of  the 
many  simplification  and  decentralized  control  methods 
currently  being  explored.  One  particular  type  of 
decentralized  control  of  interest  in  this  investigation  is 
the  decoupled  control  as  presented  by  Calico  (1).  This  type 
of  control  is  a  system  of  individual  subcontrollers  designed 
on  the  basis  of  a  subset  of  system  dynamics  with  control 
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authority  enacted  through  the  total  system  of  structural 
actuators.  Transforms  are  applied  to  the  subcontrollers  to 
insure  the  dynamic  coupling  called  observation  spillover  and 
control  spillover  is  not  destabilizing  to  the  global  system. 

In  space  a  structure  is  required  to  perform  its  mission 
for  long  durations  of  time  with  no  maintenance.  Even  if 
future  orbital  transportation  becomes  more  available,  the 
costs  associated  with  periodic  maintenance  of  active  control 
systems  will  dictate  the  need  for  highly  reliable  control 
systems.  Recently,  work  has  been  directed  towards  designing 
the  control  system  to  be  fundamental ly  more  reliable  based 
upon  its  configuration.  The  design  of  space  power  systems 
has  traditionally  used  redundant  components  in  parallel  as  a 
means  to  increase  reliability  through  system  configuration. 

A  control  system  could  in  a  similar  manner  use  redundant 
actuators,  sensors,  etc  to  build  a  more  reliable  system,  but 
it  could  also  use  more  subtle  methods.  Carignan  (4:12) 
demonstrates  that  a  control  system's  total  reliability  will 
change  based  upon  the  number  and  the  location  of  sensors  and 
actuators.  While  Carignan' s  global  searching  methods  were 
generally  applied  to  centralized  controllers  it  leads  to 
another  type  of  reliability  issue  associated  with  decoupled 
controllers.  An  individual  subcontroller  observes  a 
structure  and  enforces  control  on  the  structure  differently 
for  different  modal  arrangements.  It  is  possible  to  smartly 
arrange  the  decoupled  controller  so  that  it  is  more  robust  to 
certain  types  of  environmental  noises,  specifically  sensor 


errors.  The  basic  errors  in  sensor  measurements  will 
increase  over  time  as  the  equipment  degrades  due  to  the  harsh 
space  environment.  A  decoupled  controller  with  a  basic 
configuration,  robust  to  sensor  noise  will  enact  adequate 
control  over  the  structure  for  a  longer  period  of  time  as  the 
sensors  degrade  and  input  incorrect  information.  This  is  a 
subtle  way  of  increasing  the  overall  reliability  of  the 
system. 

The  techniques  used  to  design  decoupled  controllers 
can't  be  classified  as  rigid  or  even  as  well  defined,  but 
there  exists  some  basic  procedures  for  selecting  design 
parameters.  One  of  the  major  design  parameters  of  a 
decoupled  controller  is  the  arrangement  of  modes  between 
individual  subcontrollers.  Sumner  (11: sec  8-2)  investigated 
selecting  the  arrangement  using  modal  cost  analysis 
techniques,  modal  internal  balancing  techniques,  and 
ascending  order  ranking.  He  concluded  that  the  ascending 
order  ranking  provided  the  simplest,  yet  best  performing 
system.  He  selected  the  performance  criteria  as  time  to 
dampen  the  structure  to  a  disturbance  and  didn't  consider 
sensor  noise  robustness.  The  major  thrust  of  this 
investigation  was  to  try  and  develop  design  techniques  to 
select  parameters  of  a  decoupled  controller  so  that  the 
system  is  robust  to  sensor  noise  and  delivers  adequate 
performance.  The  main  design  parameters  investigated  were 
subcontroller  modal  arrangement  and  controller  gain. 

This  report  first  describes  the  structural  model  that 


was  used  to  demonstrate  the  design  techniques.  Next,  the 
decoupled  control  theory  is  developed.  A  chapter  is  then 
devoted  to  explaining  the  basic  statistical  tools  used.  The 
Investigation  Chapter  then  demonstates  the  actual  design 
techniques  used  by  incorporating  the  theory  developed  in  the 
previous  three  chapters.  A  Results  Chapter  then  evaluates 
the  new  controller  design  against  a  traditional  controller 
design.  The  final  chapter  summarizes  the  investigation  and 
recommends  areas  of  further  study. 
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I I .  System  Model 


All  new  structural  control  design  techniques  developed 
and  refined  in  this  thesis  were  applied  to  a  single  model 
that  simply,  but  adequately,  represents  large  space 
structures.  The  Charles  Stark  Draper  Laboratory  Inc.  of 
Cambridge  Massachusetts  test  model  one  (CSDL-I)  was  selected 
as  the  test  system  model  because  of  its  simplicity,  lightly 
damped  3-dimensional  characteristics,  and  its  general 
familiarity.  The  CSDL-I I  model  is  a  different,  more 
complicated  structure,  than  the  CSDL-I  consisting  of  59  nodes, 
numerous  truss  elements,  mirrors  and  solar  panels.  The 
CSDL-II  is  often  used  in  current  control  studies,  but  was  not 
chosen  here  because  of  its  complexity.  The  smaller  CSDL-I 
allowed  verification  of  the  fractional  experiment  techniques 
through  full  factorial  experiments.  Obviously  the  objective 
of  developing  adequate  fractional  design  techniques  becomes 
paramount  if  any  of  this  thesis  is  to  be  applied  to  actual 
structures  with  large  number  of  vibrational  modes  like  the 
CSDL-II . 

The  following  sections  will  first  describe  the  CSDL-I 
structure,  then  develop  the  necessary  equations  of  motion. 

The  CSDL-I  is  modeled  using  finite  element  techniques  and 
then  transformed  into  modal  coordinates.  The  modal  coordinate 
representation  is  then  transformed  to  a  state-space  form 
suitable  for  control  law  implementation. 


Description  of  CSDL-I 


The  CSDL-I  is  a  lumped  mass  tetrahedral  truss  shown  in 
Figure  2-1.  The  four  one  unit  masses  lumped  at  each  vertex 
of  the  tetrahedral  were  modeled  as  having  three  translational 
degrees  of  freedom.  To  enforce  only  translational  motion 
each  truss  member  was  modeled  as  a  pin  connected  member 
capable  of  carrying  only  axial  loading.  The  truss  members 
were  considered  as  massless  springs  (the  nodes  containing  all 
of  the  mass).  The  nodes  and  truss  members  are  labeled  in 
Figure  2-1. 


The  coordinates  of  each  node  are  tabulated  in  global 
coordinate  frame  s  in  Table  2-1. 


Table  2-1 

CSDL-I  Node  Coordinates 


Node 

-1 

-2 

53 

1 

0.0 

0.0 

10.165 

2 

-5.0 

-2.887 

2.0 

3 

5.0 

-2.887 

2.0 

4 

0.0 

5.774 

2.0 

5 

-6.0 

-1.155 

0.0 

6 

-4.0 

-4.619 

0.0 

7 

4.0 

-4.619 

0.0 

8 

6.0 

-1.155 

0.0 

9 

-2.0 

5.774 

0.0 

10 

2.0 

5.774 

0.0 

The  three  lower  mass 

nodes 

2,3,4  are  supported  by  two  truss 

legs  each,  which  are 

then. 

in  turn,  anchored 

.  The  nodes  5 

through  10  were  modeled  as 

perfectly  fixed. 

The  cross- 

sectional  area's  (Ac),  length  (h)  and  the  ratio  Ey(Ac)/h  for 
each  element  are  listed  in  Table  2-2.  Ey  is  defined  as 
Young's  Modulus  and  was  assumed  to  be  unity  throughout  the 
investigation  to  simplify  computations. 

Table  2-2 

CSDL-I  Element  Properties 
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The  system  was  observed  by  six  position  sensors  mounted 
on  the  six  support  legs  7  through  12.  The  sensors  measure 
axial  displacement  along  their  respective  leg  axis  in  the 
local  coordinate  frame  i.  The  system  was  controlled  by  six 
actuators  collocated  with  the  six  position  sensors.  The 
actuators  also  exerted  their  force  along  the  respective  leg 
axes.  There  was  no  requirement  to  locate  the  sensors  and 
actuators  in  such  a  configuration,  however,  the  collocated 
arrangement  had  been  used  in  previous  studies  and  allows 
easier  verification  of  computer  programs.  It  would  be  of 
benefit  to  examine  the  effects  of  differing  configurations 
since  this  would  have  a  dramatic  impact  on  performance;  this 
was  not  investigated  in  this  project  due  to  time  constraints. 
There  is  a  specific  requirement  as  to  the  number  of  sensors 
and  actuators  required  depending  on  the  particular  controller 
arrangement  selected  as  discussed  in  the  next  chapter. 

The  objective  of  the  control  system  was  to  maintain  a 
"line  of  sight"  or  LOS.  The  LOS  criteria  is  defined  as  the 
deviation  of  node  1  from  the  undisturbed  position  in  the 
global  plane  s^,  s2<  This  objective  simulates  requirements 
necessary  for  aligning  solar  panels,  antennas,  LASER  beams, 
or  many  other  uses  associated  with  orbiting  vehicle 
operations.  A  number  of  metrics  were  investigated  as 
applicable  to  represent  the  LOS  requirement  are  discussed  in 
later  chapters. 

The  following  sections  develop  the  equations  of  motion 
necessary  to  model  the  CSDL-I  and  to  enact  active  feedback 


control.  The  different  coordinate/vector  spaces  should  be 
noted  carefully  in  the  development.  The  development  starts 
in  a  local  coordinate  system  i  then  proceeds  to  a  global 
coordinate  system  s,  then  a  modal  vector  space  rj,  and  finally 
orthogonal  spaces  for  each  subcontroller's  feedback. 

Development  of  the  CSDL-I  Model  Equations 

Finite  Element  Model.  The  basic  finite  element 
vibrational  model  was  constructed  using  the  assumptions  in 
the  previous  sections  and  the  methodology  outlined  by 
Meirovitch  (7:300).  A  detailed  development  is  presented  here 
since  the  transformations  necessary  to  apply  noise  to  the 
sensors  in  local  coordinates  are  the  same  as  those  contained 
in  the  finite  element  model.  Often  the  program  NASTRAN  is 
used  to  calculate  the  equations  of  motions,  but  its  use  was 
precluded  due  to  problems  associated  with  applying  the  noise 
in  a  repetitive  computational  scheme.  The  development  starts 
with  the  system  stiffness  matrix,  then  the  mass  matrix,  then 
damping  matrix  and  finally  ends  with  the  force  input  location 
matrix . 

By  examining  a  basic  truss  element  in  a  local  coordinate 
frame,  Figure  2-2,  a  relationship  is  developed  relating  nodal 
displacement  vector  r  to  nodal  force  vector  f. 


Figure  2-2  Element  Truss  Member 


By  assuming  constant  cross-sectional  area  and  Young's  Modulus 
the  axial  displacement  has  an  equation  of  the  following  form 


0  £  x  £  h  (2.1) 


where 


Ey  =  Young’s  Modulus 

Ac  =  cross-sectional  area 

r(x)  =  axial  displacement 

h  =  length  of  truss  member 

Integrating  Eq  (2.1)  and  using  the  fact  that  the  displacement 
at  x=0  is  and  the  displacement  at  x=h  is  r2  to  solve  for 
the  constants  of  integration,  the  following  expression 
results: 


r(x)  =  (l-x/h)r^  +  (x/h)r2  (2.2) 

where 

r^  =  displacement  at  node  1 
r2  =  displacement  at  node  2 

The  displacement  equation,  Eq  (2.2),  can  be  differentiated  and 
substituted  into  the  following  boundary  conditions: 
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dr(x) 


(2.3) 


Ey(Ac) 

dx 


x=0 


dr(x) 

Ey (Ac) - 

dx 


x=h 


-2 


where 

=  force  at  node  1 
fj  =  force  at  node  2 
The  resulting  equations  are: 


(2.4) 


E2  “  £i 

Ey(Ac)  — - - —  =  -fj 


(2.5) 


Ey(  Ac )  — — - i  =  f2 


(2.6) 


Eqs  (2.5)  and  (2.6)  may  be  written  in  a  matrix  form  as: 


Kr  =  f 


where 


K 


f  1 


Ey(Ac) 


0 

0 

0 

0 

0 

0 


0-10 
0  0  0 

0  0  0 

0  10 

0  0  0 

0  0  0 


0  1 


0 

0 

0 

0 

0 


(2.7) 


(2.8) 


(2.9) 


(2.10) 
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The  stiffness  matrix  K  for  an  individual  element  can  be 
assembled  with  the  other  element  stiffness  matrices  to  form 


the  overall  system  stiffness  matrix.  The  combining  of  the 
matrices  must  be  done  in  a  global  coordinate  frame;  e.g. 
reference  frame  s  already  introduced  in  Figure  2-1.  For 
development  purposes  let  i  be  the  local  coordinate  frames  for 
each  element.  As  shown  in  Figure  2-3,  the  i^  direction 
vector  is  always  aligned  along  the  axial  direction  pointing 
from  the  lower  node  number  to  the  higher  node  number. 


Figure  2-3  Local  Coordinate  Frame 


Define  a  transform  matrix  L  such  that 


r.  =  Lr  (2.11) 

-i  -s 

where 

r^  =  the  displacement  vector  in  local  coordinate  frame  i 

r  =  the  displacement  vector  in  global  coordinate  frame  s 
”S 


L 


(transform  matrix) 
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It  should  be  observed  that  Ln=L^2  an<*  *s  t^ie  ro^ation  matrix 
of  the  global  coordinate  frame  to  the  local  frame.  That  is 

i  =  L^s  (2.12) 

Rotation  matrix  can  be  constructed  by  taking  the  product 

of  the  direction  cosines  matrices  for  each  axis  of  rotation 
necessary  to  align  the  element  axis  up  with  the  s^  direction 
vector.  By  defining  the  following  matrices  the  appropriate 
rotations  for  the  CSDL-I  can  be  tabulated  as  is  done  in  Table 
2-3: 


R2  = 


1  0 

0 

0  cos <p 

-sin# 

1  sin# 

cos# 

cosy  0 

siny 

0  1 

0 

-siny  0 

cosy 

cos# 

sin#  0 

sin# 

cos#  0 

0 

0  1 

(2.13) 


(2.14) 


(2.15) 


where 


=  rotation  angle  about  axis  i, 
=  rotation  angle  about  axis  i, 


=  rotation  angle  about  axis  i. 
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Table  2-3 

Summary  of  Rotation  Matrices 


Element 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


Rotation  Matrix 

R3(1506)R2(-54.73°) 

I 

R3(30°)R2(-54.73°) 

R-(-90°)R-(-54.73°) 

R3(-60-> 

R3(120°) 

R3(-120°)R2(-45°) 

R3(606)R2(-45°) 

R3(120°)R2(-45°) 

R,(-60“)R,(-45°) 

R3(180°)R2(-45‘) 


The  stiffness  matrix  K  and  the  force  vector  f  can  be 
converted  to  the  global  coordinate  frame  by  the  following 
transformations  developed  in  the  discussion  by  Meirovitch 
(7:315) . 


K  =  LtK.L 

S  1 


(2.16) 


-fs  =  L^i 


(2.17) 


where 


K  =  stiffness  in  global  coordinate  frame  s 
s  - 

=  stiffness  in  local  coordinate  frame  i 
f  =  force  in  global  coordinate  frame  s 

“  S  “ 

f ^  =  force  in  local  coordinate  frame  i 
After  transforming  each  of  the  12  stiffness  matrices 
associated  with  the  12  truss  elements  to  the  global 
coordinate  frame  s  they  are  summed  together  in  their  extended 
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form.  The  extended  form  is  a  30x30  matrix  where  zero  rows 
and  zero  columns  have  been  inserted  in  the  original  6x6  K 

8 

matrices  for  those  generalized  nodal  displacements  that  are 

not  represented  in  the  original  K  .  An  example  extended 

s 

matrix  for  element  3  is  listed  in  Appendix  A.  After  summing 

the  extended  K  matrices,  the  boundary  conditions  are  applied, 
s 

The  boundary  conditions  are  enforced  by  deleting  the  rows  and 
columns  associated  with  the  fixed  nodal  displacements  (i.e. 
nodes  5-10)  since  reaction  forces  are  not  of  interest.  The 
final  stiffness  matrix  associated  with  the  equation  of  motion 
is  also  included  in  Appendix  A.  The  remainder  of  this  report 
shall  refer  to  the  total  12x12  stiffness  as  simply  K  unless 
noted  otherwise.  A  computer  program  to  calculate  K  called 
STIFF. M  written  in  Pro-MATLAB  script  language  is  listed  in 
Appendix  B. 

The  mass  matrix  M  is  simply  an  identity  matrix  I  since 
the  original  assumption  was  1  unit  lumped  masses  at  each  of 
the  nodes.  It  should  be  noted  that  if  the  assumption  of 
lumped  masses  was  not  enforced  the  mass  matrix  should  be 
developed.  The  mass  matrix  can  be  constructed  as  outlined  in 
Meirovitch  (7:306)  using  the  same  interpolation  function 
(i.e.  Eq  (2.2))  as  the  stiffness  matrix.  The  method  is  not 
discussed  here  since  it  was  not  used. 

The  general  form  of  the  system  equation  of  motion  may  be 
written  in  the  generalized  global  coordinate  frame  as: 


where 


M  =  nxn  mass  matrix 
E  =  nxn  damping  matrix 
K  =  nxn  stiffness  matrix 

D  =  nxm  input  location  matrix  (actuator  locations) 
s  =  nxl  global  coordinate  vector 
u  =  nxl  control  input  vector 

Matrix  E  contains  all  of  the  viscous  damping  terms.  The 
CSDL-I  model  does  not  incorporate  any  assumptions  about  the 
damping.  The  viscous  terms  could  be  modeled  in  the  global 
coordinate  frame  and  transformed  into  modal  coordinates,  much 
like  the  stiffness  terms.  The  disadvantage  is  that  the  E 
matrix  will  not  be  decoupled;  it  is  more  convenient  to  assume 
a  representative  viscous  damping  term  which  can  be  modeled 
directly  in  modal  coordinates.  An  assumed  value  of  .005  for 
each  mode  was  use  to  represent  a  lightly  damped  system. 

The  input  location  matrix  D  converts  the  force  applied 
from  any  one  actuator  to  the  global  coordinate  frame.  It 
should  be  noted  that  u  isn't  in  the  global  coordinate  frame, 
but  represents  a  vector  of  actuator  gains  along  their  own 
local  coordinate  axis.  The  D  matrix  is  constructed  by 
transforming  the  actuator's  i^  direction  vector  to  global 
coordinate  frame  using  the  relationships  in  Eqs  (2.12)  and 
(2.17).  The  matrix  is  then  assembled  like  the  K  matrix  using 
an  extended  form.  A  program  to  calculate  D  called  SENSOR. M 
written  in  Pro-MATLAB  script  language  is  listed  in  Appendix 
C.  This  code  is  used  extensively  to  calculate  the  final 
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state-space  matrices  after  applying  perturbations  to  the 
actuators  in  the  local  coordinate  frame.  The  final  D  matrix 
is  listed  in  Appendix  A. 

Modal  Coordinate  Transformation.  The  modal  coordinate 
system  is  a  vector  space  which  utilizes  the  natural  mode 
shapes  as  a  basis.  These  natural  modes  are  orthogonal  to 
each  other  and  therefore  will  decouple  the  equations  of 
motion.  This  decoupling  will  allow  each  individual  second 
order  differential  equation  to  be  represented  as  two  coupled 
first  order  differential  equations.  The  two  coupled  first 
order  differential  equations  allow  for  conversion  to 
state-space  form,  useful  for  applying  structural  control 
techniques . 

The  natural  mode  shapes  and  their  associated  frequencies 
can  be  computed  by  solving  the  characteristic  eigenvalue 
problem. 

[K  -  XjMHUj}  =  0  (2.19) 

where 

Xj  *  eigenvalue 

Uj  =  eigenvector 

It  can  also  be  shown  that  the  natural  frequencies  of  the 
system  are: 


where 

2 

=  natural  frequency 

A  complete  description  is  provided  in  Mierovitch  (7:107). 


Table  2-4  lists  the  natural  frequencies  for  the  CSDL-I  model. 


Table  2-4 

CSDL-I  Natural  Frequencies 

Assigned  Mode  Natural  Frequency 

Number _ o  in  Rad/sec 


1 

1.34 

2 

1.66 

3 

2.89 

4 

2.96 

5 

3.40 

6 

4.20 

7 

4.66 

8 

4.75 

9 

8.54 

10 

9.25 

11 

10.28 

12 

12.90 

Once  the  eigenvalue  problem  is  solved,  U  becomes  a  matrix 
whose  columns  are  the  IK  eigenvectors  associated  with  the  Xj 
eigenvalues.  It  should  be  noted  that  the  amplitudes  of  IK 
are  not  unique,  only  the  shape.  A  unique  amplitude  can  be 
enforced  by  requiring  orthonormality  of  U  with  the  mass 
matrix  M. 

UtMU  =  I  (2.21) 

The  equations  of  motion  can  then  be  transformed  to  modal 
coordinates  through  the  relationship 

s  =  Urj  (2.22) 

where 

rj  =  the  vector  of  modal  coordinates 
Substituting  Eq  (2.22)  into  Eq  (2.18)  yields 

UTMrj  +  UTEUn  ♦  UTKrj  =  UTDu  (2.23) 

By  first  recognizing  U  is  a  matrix  of  orthogonal 


eigenvectors,  M  is  an  identity  matrix,  and  using  the 
relationship  in  Eq  (2.19),  Eq  (2.23)  can  be  rewritten  as: 


if}  +  UTEUrj  + 


UtDu 


(2.24) 


Recalling  that  the  viscous  damping  is  modeled  only  in  the 
modal  coordinates,  Eq  (2.24)  can  be  expressed  as: 


Jri  + 

^  2 
2C<0j 

D  + 

N  2 

Wj 

\ 

\ 

rj  =  U  Du 


(2.25) 


where 


£  =  modal  damping  (assumed  as  .005) 

State-space  Transformation.  After  obtaining  the 
equations  of  motion  in  modal  coordinates,  Eq  (2.25),  a 
transformation  to  a  state  space  representation  is  possible. 
As  presented  in  Sumner  (11: Sec  3,1)  the  objective  is  to 
arrange  the  equations  into  the  form: 


x  =  Ax  +  Bu  (2.26) 

y  =  Cx  (2.27) 


where 

A  =  2nx2n  plant  matrix 
B  =  2nxm  control  input  matrix 
C  =  2nxl  output  matrix 
x  -  2nxl  state  vector 
u  =  mxl  control  input  vector 

This  conversion  is  implemented  by  defining  the  state  vector 


This  in  turn  requires  the  other  matrices  to  be  of  the  form 


(2.29) 


(2.30) 


c  =  tcpu  S  CyU]  (2.31) 

where 

C  =  n  xn  displacement  sensor  output  matrix 
P  P 

Cv  =  nyxn  velocity  sensor  output  matrix 
In  the  particular  arrangement  of  the  CSDL-I  model  used  in 
this  report  there  are  no  velocity  sensors,  so  Cv  is  simply  a 
matrix  of  zeros.  Also,  the  sensors  are  collocated  with 
actuators  so  that  the  Cp  matrix  is  simply  the  transpose  of 
the  D  matrix. 

C  =  DtU  or  [UtD]t  (2.32) 

P 

The  final  output  matrix  for  the  CSDL-I  is  then 

C  =  C(UtD)t  |  0  ]  (2.33) 

The  relationship  in  Eq  (2.33)  is  used  to  invoked  noise  in  the 
sensors  during  feedback.  As  shown  in  later  sections 
perturbations  are  applied  in  local  i  coordinates,  then  a  new 
D  matrix  is  calculated,  and  then  finally  a  new  C  matrix  is 


formed  to  be  used  in  the  control  methodol 
B,  and  C  matrices  are  listed  in  Appendix 
develops  the  control  methodology  applied 
Eqs  (2. 28), (2. 29), (2. 30), and  (2.33). 


ogy.  The  nominal  A, 
D.  The  next  chapter 
to  the  state-space 
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Ill .  Control  Law  Development 


Decoupled  control  can  be  implemented  with  a  variety  of 
different  dynamic  feedback  control  laws.  The  decoupled 
methodology  uses  orthogonal  transforms  to  decompose  the 
original  single  state-space  representation  into  a  number  of 
separate  representations.  These  separate  state-space 
representation  can  have  their  own  individual  dynamic  feedback 
control  laws  applied  to  them  without  destabilizing  the  other 
individual  systems.  It  should  be  noted  that  any  unmodeled 
dynamics  (i.e.  modes)  may  or  may  not  be  destabilized  by 
applying  dynamic  feedback  control  to  the  separate 
subcontrollers.  This  fact  is  true  for  any  dynamic  feedback 
system,  decoupled  or  not,  that  has  only  partial  knowledge  of 
the  plant  dynamics  within  the  bandwidth  of  the  control 
system. 

The  first  of  the  following  sections  introduces  the 
fundamentals  of  linear  quadratic  feedback  theory.  This 
particular  type  of  dynamic  feedback  was  implemented  on  a 
complete  fullstate  representation  of  the  CSDL-I  as  a 
benchmark  for  decoupled  controller  performance.  The  linear 
quadratic  regulator(LQR)  and  linear  quadratic  estimator(LQE) 
theories  were  also  applied  to  the  separate  decoupled 
state-space  systems.  The  LQR/LQE  theory  is  developed 
separately  to  emphasize  that  other  types  of  dynamic  feedback 
control  laws  could  be  implemented.  The  second  section 
outlines  the  theory  necessary  to  decouple  the  original 


state-space  representation  into  separate  decoupled  systems. 

The  third  section  then  applies  the  theory  of  the  first  two 
sections  to  the  CSDL-I. 

Full state  Linear  Quadratic  Theory 

As  presented  more  thoroughly  in  Franklin  (6:312-405), 
Stunner  (ll:sec  3,6-9)  and  Ridgley  (9:sec  6,1-15),  a  state-space 
system  of  the  form 

x  =  Ax  +  Bu  (3.1) 

y  =  Cx  (3.2) 

may  be  actively  controlled  by  "feeding  back”  a  set  of  control 
forces  based  upon  the  current  state  of  the  system,  or 
mathematically  let 

u  =  Kcx  (3.3) 

where 

Kc  =  matrix  of  control  gains 

Eq  (3.3)  assumes  complete  knowledge  of  the  state  vector  x. 

Since  this  is  not  the  case  with  the  given  six  position 
sensors  (x  is  rarely  known  completely),  a  method  to  estimate  x 
is  constructed.  The  x  vector  can  be  estimated  by  taking 
information  from  the  sensors  y  and  forming  an  "estimator"  or 
sometimes  called  "observer"  of  the  form 

ft  =  Aft  +  Bu  +  Kf {y-ft}  (3.4) 

ft  =  eft  (3.5) 

where 


£  =  estimated  state  vector 

ft  =  estimated  output  vector 

Kf  =  matrix  of  estimator  gains 

The  feedback  is  then  accomplished  using  the  estimated  states. 

u  =  Kcft  (3.6) 

The  Kc  and  Kf  matrices  are  formed  using  LQR/LQE  theory.  The 
Kf  matrix  must  also  be  of  a  form  to  insure  the  error  state 
defined  as: 

e  =  ft  -  x  (3.7) 

where 

e  =  state  error  vector 
goes  asymptotically  to  zero  in  steady  state. 

The  LQR  theory  determines  Kc  by  minimizing  a  cost 
function  J  that  represents  a  compromise  between  obtaining 
minimum  deviation  of  the  present  state  from  the  final  desired 
value  (tracking  error)  and  control  input  effort  necessary  to 
obtain  the  desired  value. 

CO 

J  =  J  (xTQ  x  +  uTRQu)dt  (3.8) 

Jo 

where 

Qo  =  a  positive  semidefinite  square  state  weighting 
matrix 

R  =  a  positive  definite  square  control  weighting  matrix 
o 

The  solution  to  the  minimization  of  Eq  (3.8)  can  be  found  to 


where 


P  =  positive  definite  solution  to  the  following  Riccati 
Equation: 

ATP  +  PA  -  PBR^1BTP  +  Qq  =  0  (3.10) 

Using  the  principle  of  duality,  the  gain  matrix  Kf  is  also 

obtained  by  the  solution  method  outlined  in  Eqs  (3.8)  (3.9) 

and  (3.10)  by  substituting  CT  for  B  and  AT  for  A.  The 

weighting  matrices  Q  and  R  are  not  necessarily  the  same  for 

o  o 

the  Kf  calculation.  The  Q  and  R  used  in  the  estimator 

o  o 

calculations  contain  assumptions  about  how  the  process  noise 
and  measurement  noise  are  distributed  across  the  state-space 
model.  To  insure  stable  feedback  for  fullstate  control 
matrices  A  and  B  must  be  a  controllable  pair  along  with  the 
conditions  on  Qq  and  Rq  stated  in  Eq  (3.8).  A  more  detailed 
discussion  is  given  in  Ridgely  (9:sec  6,5-8).  This 
investigation  used  diagonal  matrices  for  Qq  and  Rq,  thus 
insuring  a  stable  LQR/LQE  compensator. 

Now  that  the  components  have  been  developed  the  complete 
compensator  design  can  be  formed.  By  substituting  Eq  (3.6) 
into  Eq  (3.1) : 

x  =  Ax  +  BKcft  (3.11) 

Substituting  Eq  (3.2),  (3.5)  and  (3.6)  into  Eq  (3.4)  yields: 

ft  =  Aft  +  BKcft  +  Kf (Cx-cft)  (3.12) 

Differentiating  Eq  (3.7)  yields: 


Substituting  Eq  (3.11),  (3.12)  into  (3.13)  collecting  terms 
and  using  Eq  (3.7)  to  define  e  yields: 

e  =  (A-KfC)e  (3.14) 

Taking  Eq  (3.7)  and  reforming  as: 

$  =  e  +  x  (3.15) 

and  substituting  into  Eq  (3.11)  results  in: 

x  =  ( A+BKc ) x  +  BKce  (3.16) 

Now  Eqs  (3.16)  and  (3.14)  can  be  written  in  matrix  form  as 


x 

e 


A  +  BKc  BKc 

0  A  -  KfC 


X 

e 


(3.17) 


or  more  compactly 


x 

e 


=  Ae 


x 

e 


(3.18) 


where 

Ae  =  the  state-space  closed  loop  plant  matrix  with  error 

states 

The  compensator,  Eq  (3.17),  doesn't  have  an  error  state 
reference  input  matrix.  Be.  Franklin  (6:372)  develops 
various  methods  of  forming  the  error  state  input  reference 
matrix.  Be,  based  on  whether  the  error  states  are  to  be 
excited  or  not.  This  investigation  looked  only  at  a 
reference  impulse  as  the  test  input,  where  an  impulse  can  be 
thought  of  as  a  set  of  initial  displacements.  The  original  B 
matrix  can  then  be  used  to  convert  the  set  of  initial 


displacements  along  the  actuator's  local  coordinate  frame  i 

•  into  the  state-space  frame,  by  defining  an  error  state  input 
matrix  in  the  form  shown  in  Eq  (3.19).  However,  using  the 
original  B  matrix  to  form  the  Be  matrix  does  restrict 

•  applying  the  initial  displacement  conditions  to  combinations 
of  displacements  along  the  actuator  axes,  but  this  presented 
little  problem  in  the  investigation.  A  typographical  mistake 

•  in  the  program  DECOUPLE. M,  Appendix  P,  was  noticed  late  in 
this  investigation  and  prevented  implementation  of  Be  as 
previously  defined.  Further  comments  on  the  necessary 

•  changes  required  in  the  code  to  implement  the  defined  Be  are 
included  in  the  program  documentation.  For  the  purposes  of 
this  investigation  the  exact  form  of  Be  was  of  little 

•  consequence  since  the  reference  inputs  used  were  exactly  the 
same  for  all  design  configurations  considered.  The  main 
objective  was  to  compare  different  design's  relative 

•  performance,  but  not  to  predict  actual  performance  of  a  real 
structure  to  a  specific  input.  The  error  state  output  matrix 
Ce  is  simply  the  original  C  matrix  augmented  with  zeros  as 

•  shown  below  since  the  sensor  locations  have  not  changed.  The 
final  compensator  design  with  LQR/LQE  control  can  be 
represented  as: 


(3.19) 
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y  =  Ce 


x 

e 


(3.20) 


where 


Ae  = 


Be  = 


Ce  = 


defined  in  Eq  (3.17) 
B 
0 

[C  •  0  ] 


Decoupled  Controller  Theory 

The  basic  concept  of  decoupled  control  is  to  break  an 
entire  structural  model  Eqs  (3.1)  and  (3.2)  into  partitions  that 
can  be  actively  controlled  separately  without  destabilizing 
any  other  partition.  This  decoupled  controller  theory 
section  also  addresses  a  partition  of  residual  structural  modes 
that  are  not  considered  in  the  control  development,  but  are 
modeled  in  the  time  response  of  controller  performance. 

These  residuals  mathematically  represent  the  partition  of  the 
physical  model  that  is  unknown  for  control  law  development, 
yet  ever  so  present  in  application. 

As  presented  in  Calico  (2:10-11)  the  state-space  vector 
can  be  partitioned  into  separate  subcontroller  and  a  residual 
components . 


x  = 


-1 

?2 


x 

-r 


(3.21) 
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where 

xN  =  a  vector  of  modes  grouped  into  the  Nth 
subcontroller 

;r  =  the  residual  modes  that  will  not  be  actively 
controlled 

The  decision  of  which  structural  modes  to  group  into  which 
subcontroller  group  x^  is  one  of  the  main  thrusts  of  this 
investigation.  It  will  be  shown  later  the  decoupling 
transform  enforces  different  restrictions  on  subcontroller  x^ 
than  xR  ,  so  that  the  performance  of  the  total  system  will  be 
different  for  different  arrangements.  Each  subcontroller  can 
be  written  as  a  partition  of  the  total  plant. 

x  =  A^x  +  B^u  i=l,...N  (3.22) 

and 

5r  =  Ar?  +  Br¥  (3.23) 

where  A^  ,  B^  ,  Ar  and  Br  are  defined  as  the  partitions  of  Eqs 
(2.29)  (2.30)  and  (2.31). 


Ai  = 


0  j  I 

\  o  i~\ . 2 . 

\ 


(-o.2). 

v  "-j  /x 


(3.24) 


Bi  = 


(UTD.). 


(3.25) 


Ct  *  [(C  D)i  1  (CyU).] 


(3.26) 
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where 

j  =  the  individual  mode  assigned  to  the  ith  controller 
i  =  the  ith  partition  of  residual  partition 
The  dimensions  of  the  partitioned  matrices  are  2n^x2n^  for  A^  , 
2n.xn  for  B.  ,and  n  x2n.  for  C.  where  n.  is  the  number  of 

1  a  1  SI  X  X 

modes  in  the  ith  subcontroller,  n  is  the  total  number  of 

8 

actuators  for  the  structure,  and  n  is  the  number  of  sensors 

s 

for  the  structure. 

To  apply  control  to  the  N  individual  systems  defined  as: 

?i  =  Aixi  +  Biu  (3.27) 

y  =  Cx  (3.28) 


requires  that  the  controllers  don't  interact  with  each  other. 
In  the  form  of  Eqs  (3.27)  and  (3.28)  interaction  is  present 
in  that  the  control  input  u  and  the  output  vector  y  are  the 
same  for  all  N  subcontrollers.  This  interaction  is  called 
control  spillover  and  observation  spillover  respectively  and 
may  be  eliminated  by  enforcing  certain  constraints  on  the  Kc 
and  Kf  matrices.  Consider  a  new  control  variable  and 

output  variable  w^  defined  as: 

N 

>  S  I  Ti-i  (3.29) 

i  =  l 


^i 


(3.30) 


By  substituting  Eq  (3.29)  into  Eq  (3.27)  and  Eq  (3.28)  into 
Eq  (3.30): 


N 


x. 

-i 


Ai?i 


+  B 


j=l 


T  j”  j 


(3.31) 
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(3.32) 


N 

si  =  ri  2  Vi 

j=i 

To  completely  remove  all  spillover  effects  and  I\  form  the 
null  space  of  the  other  separate  partitions,  and  matrices, 
or  symbolically. 


=  0 

i=l .  .  .N 

for 

j=l. .  .N 

(3.33) 

riCj  '  0 

It  is  also  possible  to  enforce  less  restrictive  conditions  on 
Ti  and  I\that  allow  coupled  but  stable  subcontrollers.  They 
are: 


B.T .  =  0  i=l . . .N-l 

1  3  for  i#j  (3.34) 

=  0  j=l+l.  .  .N 

The  conditions  in  Eq  (3.34)  were  used  in  this  investigation 
to  enforce  a  stable  set  of  subcontrollers.  Now  noticing  that 
by  substituting  conditions  in  Eq  (3.34)  into  Eqs  (3.31)  and 
(3.32)  the  following  maybe  formed. 

x.  =  Ai?i  +  Bi-i  (3.35) 


Si  *  Ci?i 


(3.36) 


where 


Bi  8  BiTi 
Ci  =  riCi 

Any  type  of  dynamic  feedback  control  can  be  implemented  on  the 
N  set  of  equations  Eq  (3.35)  and  Eq  (3.36)  without 
destabilizing  each  other.  The  LQR/LQE  theory  as  presented  in 
the  previous  section  was  implemented  on  Eqs  (3.35)  and  (3.36) 
resulting  in  N  controller  gain  matrices  Kc  and  N  estimator 
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gain  matrices  Kf.  The  notation  of  Kc^  and  Kf^  will  be  used 
to  denote  controller  and  estimator  gain  matrices  for  the  ith 
individual  subcontroller.  The  resulting  form  of  the 
estimator  for  each  subcontroller  as  defined  in  Eqs  (3.35)  and 
(3.36)  is: 

i.  =  A.$  +  B*$.  +  Kf.(w.-$.)  (3.37) 

where 

«i  -  <Si 

and  for  feedback  control 

9±  =  Kcifti  (3.38) 

To  implement  the  subcontrollers,  the  compensator  system  as 
defined  in  Eqs  (3.37)  and  (3.38)  are  transformed  back  to  a 
form  containing  the  actual  input  vector  u.  The  vector  u 
defines  the  set  of  actuators  that  may  physically  be  invoked 
to  perform  control.  Substituting  Eq  (3.38)  into  (3.39)  yields 
feedback  control  equations: 

Ji  =  A.x.  +  B^u  (3.39) 

N 

u  =  J  TiKci$i  (3.40) 

i  =  l 

The  above  is  another  expression  of  Eq  (3.31)  that  indicates 
that  u  will  be  the  input  control  vector  instead  of  $.  A  new 
closed  loop  plant  matrix  Ae  can  be  formed  by  applying  the 
steps  outlined  for  a  fullstate  controller  Eqs  (3.11)  through 
(3.17)  to  each  of  the  N  subcontrollers  and  then  augmenting 
the  individual  Ae^  into  a  single  Ae  that  also  includes 


3-11 


residuals.  Instead  of  developing  the  general  form  which  » 

slightly  messy  the  procedure  will  be  applied  directly  to  the 

CSDL-I  model  in  the  the  following  section. 

Before  developing  the  actual  model  the  number  of  sensors 

and  actuators  required  to  enforce  decoupled  control  are 

examined.  The  conditions  in  Eq  (3.34)  to  enforce  orthogonal 

subcontrollers  require  the  column  of  to  be  orthogonal  to 

the  rows  of  B1  through  BN_^.  Calico  (2:11)  shows  that  for 

to  exist  there  must  be  enough  actuators  to  make  B^  be  of 

sufficient  rank  to  ensure  a  null  space  of  at  least  one  vector 

exists.  The  exact  requirement  is 
N-l 

na  >  l  ni  (3'41) 

i  =  l 

In  words  this  says  the  number  of  modes  grouped  in  the  first 

N-l  subcontrollers  must  be  less  than  the  available  actuators. 

Similarly  in  order  for  to  exist  there  must  be  enough 

sensors  to  make  Ci  be  of  sufficient  rank.  That  is: 

N-l 

ns  >  ^  nA  (3.42) 

i=l 

Again,  in  words,  this  states  the  number  of  sensors  must  exceed 
the  number  of  modes  grouped  in  the  last  N-l  subcontrollers. 

Controller  Model  for  CSDL-I 

The  particular  arrangement  of  subcontrollers  used  in 
this  investigation  involved  three  modes  grouped  in 
subcontroller  one,  two  modes  grouped  into  subcontroller  two, 
and  three  modes  grouped  into  subcontroller  three,  with  the 
highest  four  modes  always  assigned  to  the  residuals.  With  an 


arrangement  of  six  sensors  and  six  actuators  the  requirements 
of  Eqs  (3.41)  and  (3.42)  are  satisfied.  The  variation  in 
size  of  the  groupings  was  not  investigated,  only  the 
assignment  of  modes  1-8  to  the  3-2-3  subcontroller 
arrangement.  It  would  be  desirable  to  limit  the  number  of 
subcontrollers  in  an  actual  implementation  so  as  to  minimize 
energy  input  necessary  to  keep  subcontrollers  orthogonal. 
Therefore,  a  realistic  design  would  choose  subcontrollers  as 
large  as  possible  based  on  computational  capabilities,  thus 
making  the  size  variation  an  unimportant  variable.  This 
decision  was  part  of  the  original  outscoping  and  downscoping 
of  the  problem. 

The  x  vector  is  first  partitioned  into  its  components 
and  augmented  with  the  error  states  to  form  the  augmented 
state  vector  z. 


*  -  »T  *T  „T  -T  T  T  ,1 

5  “  •  SjL  *  ?2  9  - 2  9  —  3  *  S3  9  5r  / 


(3.43) 


The  subcontroller's  states  and  feedback  are  formed  using  Eqs 
(3.39)  and  (3.40). 


-1  =  Al-1  +  Bl- 


(3.44) 


x2  =  AjXj  +  BjU 


(3.45) 


x3  =  A3x3  +  B3u 


(3.46) 


?r  =  Ar-r  +  Br- 


(3.47) 


U  =  Tj^KCj^  +  TjKCj  +  TgKCgl, 


(3.48) 
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The  estimators  are  formed  by  applying  Eq  (3.37)  written  in 


input  vector  u  notation,  and  also  Eq  (3.32). 
-1  =  Al$l  +  Bl-  +  Kfl*-1~-1* 


(3.49) 


X,  =  A,$,  +  B2u  +  Kf2(w2~02) 


(3.50) 


=  A,$,  +  B3u  +  Kf3(w3~03) 


(3.51) 


-1  =  riCl?l  +  riC2-2  +  riC3-3  +  riCr5r 


(3.52) 


-2  =  F2C1-1  +  r2C2?2  +  r2C3?3  +  r2Cr5r 


(3.53) 


-3  =  r3Cl?l  +  r3C2-2  +  r3C3?3  +  r3Cr?r 


(3.54) 


$1  ■ 


(3.55) 


-2  =  r2C2-2 


(3.56) 


-3  «  r3C3£3  (3.57) 

By  substituting  Eq  (3.48)  into  Eq  (3.44)  and  Eqs  (3.52)  and 
(3.55)  into  Eq  (3.49)  and  then  recalling  ^  yields 

the  error  state  equations: 


§i  =  (Ax  -  +  Kf^C 2*2  + 


+  Kfir1Crxr 


(3.58) 


e2  =  (A2  -  Kf2r2C2)e2  +  Kf2r2Cl?l  +  Kf2F2C3?3 


+  Kf2r2Crxr 


(3.59) 


e3  =  (A3  -  Kf3r3C3)e3  +  +  Kf^CjXj 


+  Kf3r3Crxr 


(3.60) 


3-14 


Substituting  $3  =  e^  +  x.  into  Eq  (3.48)  and  then  combining 
Eq  (3.48)  with  Eqs  (3.44),  (3.45),  (3.46)  and  (3.47)  yields 
the  state  variable  portion  of  z: 


Xi  *  (A3  +  B1T1Kc1);1  +  BiTiKc1Sl  +  biT2Kc2^-2  +  -2^ 

+  8^3X03(03  +  x3)  (3.61) 

x2  =  (A,  +  B2T2Kc2)x2  +  B2T2Kc2e2  +  62X3X03(53  +  X3) 

+  B2T3Kc3(e3  +  x3)  (3.62) 

*3  =  (A3  +  +  B3T3KC3?3  +  B3TlAcl(  ®1  +  ?1> 

+  B3T2K c2(e2  +  x2)  (3.63) 

xr  =  Arxr  +  8^3X03(53  +  X3)  +  BrT2Xc2(e2  +  x2) 

+  BrT3Xc3(e3  +  ;3)  (3.64) 


Now  Eqs  (3.58),  (3.59),  (3.60),  (3.61),  (3.62),  (3.63),  and 
(3.64)  may  be  assembled  into  the  matrix  form  of  Ae,the 
compensator  closed  loop  plant  matrix: 


z  =  Aez  (3.65) 

where 
Ae  = 


Al+BlKcl 

B3XC3 

BlT2Kc2 

B1T2Kc2 

b1T3Kc3 

B1T3Kc3 

0 

0 

Al~KflCl 

Xfiric2 

0 

KflFlC3 

0 

Xf3r3Cr 

B2T1KC1 

B2T1Kc1 

* 

A2+B2Xc2 

B2Kc2 

B2T3KC3 

B2T3Kc3 

0 

Xf2r2Ci 

0 

0 

A2'Kf2C2 

Xf2r2C3 

0 

xf2r2cr 

B3T1KC1 

B3T1KC1 

B3T1KC2 

B3T1KC2 

A3+B3KC3 

B3KC3 

0 

Kf3F3Cl 

0 

Kf3r3C2 

0 

0 

Aj-Xf 3C3 

Xf3r3cr 

8^3X03 

8^3X03 

BrT2Xc2 

BrT2Xc2 

BrT3Xc3 

BrT3Xc3 

Ar 
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It  can  easily  be  seen  that  by  applying  the  requirements  of 

Eq(3.34)  that  the  Ae  matrix  becomes  lower  block  triangular. 

This  condition  suppresses  the  spillover.  An  example  Ae 

matrix  for  the  CSDL-I  is  contained  in  Appendix  E  for  modes 

1,2,3  assigned  to  subcontroller  one,  modes  4,5  assigned  to 

subcontroller  two,  and  modes  6,7,8  assigned  to  subcontroller 

three.  The  controller  and  estimator  gains  where  computed  as 

previously  explained  both  using  a  diagonal  Q  matrix  of  20 

o 

and  a  unity  matrix  R  . 

o 

The  next  chapter  will  depart  from  the  equations  of 
motion  and  develop  the  statistical  techniques  necessary  to 
efficiently  exercise  these  equations  of  motion.  These 
techniques  will  be  used  to  create  a  numerical  experiment 
plan.  The  numerical  experiment  consisting  of  calculating  the 
time  history  of  the  system  to  various  applications  of 
impulses  to  Eqs  (3.19)  and  (3.20)  using  the  definition  of  Ae 
as  given  in  Eqs  (3.17)  and  (3.65). 


IV .  Response  Surface  Methodology 


Box  and  Draper  (1:1)  explain,  "Response  surface 
methodology  comprises  a  group  of  statistical  techniques  for 
empirical  model  building  and  model  exploitation.  By  careful 
design  and  analysis  of  experiments,  it  seeks  to  relate  a 
response,  or  output  variable  to  the  levels  of  a  number  of 
predictors,  or  input  variables,  that  effect  it."  Usually 
response  surface  methodology  (RSM)  is  applied  to  experimental 
data  or  simulation  models,  but  can  be  applied  to  a 
deterministic  type  model  such  as  the  mathematical 
representation  of  the  CSDL-I  previously  developed.  The 
portions  of  RSM  used  in  this  investigation  were  principally 
regression  techniques,  fractional  designs  and  signal -to-noise 
concepts.  These  tools  allowed  intelligent  probing  of  the 
large  number  of  design  possibilities  and  identification  of 
important  design  variables.  They  three  tools,  regression 
techniques,  fractional  designs  and  signal-to-noise  concepts, 
will  be  developed  in  the  following  separate  sections  in  the 
order  mentioned. 

Linear  Regression 

By  investigating  a  complicated  relationship  over  a 
sufficiently  restricted  region  any  function  can  be 
approximated  by  a  polynomial.  This  approach  involves  the 
Taylor's  series  expansion  of  the  true  function  over  the  region 
of  interest.  Let  the  proposed  polynomial  approximation  be 
represented  as: 
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where 


y  =  output  (i.e.  response  time,  etc) 
f  =  polynomial  function 
5  =  set  of  input  variables 

0  =  set  of  parameters  or  coefficients  in  equation 
e  =  experimental  error 

The  method  of  least-squares  estimates  0  by  choosing  a  value 
of  0  such  that  the  sum  of  squares  is  minimized.  The  sum  of 
squares  S(0)  being 

s(e)  =  J  {yu-fUU'0)>2  <4*2> 

where 

yu  =  actual  output  value  at  design  point  u 

=  particular  setting  of  inputs  for  design  point  u 
The  least -squares  estimate  depends  upon  experimental  errors 
being  statically  independent,  constant  variance  and  normally 
distributed.  This  is  often  verified  by  examining  a  plot  of 
the  residuals. 

Eq  (4.1)  can  be  rewritten  in  a  matrix  form  assuming  a 
linear  function  f. 

y  =  Z0  +  £  (4.3) 


where 


nP  =  number  of  data  points 
p  =  number  of  parameters  in  assumed  model 
Z  =  experimental  design  matrix 
Using  the  form  of  Eq  (4.3)  it  can  be  shown.  Box  and  Draper 


(1:40-42),  that  the  estimate  of  0  that  minimizes  S(0)  is 


(ZTZ)-1ZTy 


(4.4) 


where 

A 

0  =  estimate  of  0 

To  investigate  a  relationship  using  least-squares  the  form  of 
/  must  be  hypothesized.  That  is  to  say  variables 
must  initially  be  considered  important.  Each  data  point  yn^ 
is  then  generated  from  the  physical  or  mathematical  model  for 
a  unique  combination  of  variables  £^...£p.  This  unique 
combination  is  called  a  design  point  (i.e.  a  row  of  matrix 
Z).  If  the  model  is  exercised  twice(replication)  at  one 
design  point  another  data  point  is  generated,  but  the  number 
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of  design  points  remains  constant.  It  is  important  to  note 
the  difference  between  design  point  and  data  point.  A 
regression  is  then  performed  on  the  data  using  Eq  (4.4).  The 
statistical  significance  of  each  8p  associated  with  a 
particular  variable  £p  and  the  overall  model  can  be  computed 
using  an  F  statistic.  This  allows  unimportant  variables  to 
be  removed  creating  a  parsimonious  model  capturing 
significant  effects  and  interactions.  The  F  statistic  is 
defined  as: 


pp  ■ 


(4.5) 


where 


F  =  F  statistic  for  parameter 
P  P 

MS  =  mean  sum  of  squares  for  parameter 
P  P 

MS  =  mean  sum  of  squares  for  pure  error  (only  possible 
e 

if  a  replication  of  a  design  point  exists) 

MS^  =  mean  sum  of  squares  for  residuals 
(lack-of-fit  when  MSe  is  zero) 


i*  MS, 
F.  =  _ t 

*  MS 

e 


(4.6) 


where 


F^  =  F  statistic  indicating  appropiateness  of  model 

MS^  =  mean  sum  of  squares  of  lack-of-fit  of  model 

The  calculations  of  MS_  and  MS  are  simplified  greatly  if  the 

®  P 

levels  of  each  £  are  chosen  so  that  the  columns  of  Z  are 

P 

orthogonal  to  each  other  (each  column  of  Z  is  associated  with 
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a  particular  variable  £p) .  This  is  discussed  in  the  next 
section.  Assuming  an  orthogonal  Z  matrix  the  following  can 
be  shown.  Box  and  Draper  (1:34-62). 

MSp  =  <ZpV>2/nP  <4*7> 

where 

Z  =  column  of  Z  associated  with  variable  £ 

p  'p 

and 


mP  rP 

MS  -ll  <Viu-Vi)2  (4.8) 

*  i=l  usl 

(nP  -  mP )  ~ 

where 

mP  =  number  of  design  points  (i.e.  unique  combinations  of 
variables  ) 

rP  =  number  of  replications  at  design  point  i 

g^  =  mean  across  all  replications  at  design  point  i 

and 


ms4  = 


■  p 

l  (zTy)2/nP 

-i  =  l _  • 

(mP  -  p) 


-2 


-  npg  -  MS  (np-mp) 


(4.9) 


where 

g  =  mean  of  all  data  points 

By  noticing  that  if  there  are  no  replications  of  design 
points  (rp  =  0)  Eq  (4.9)  can  be  simplified  because  MS 

© 

reduces  to  zero. 


4-5 


(4.10) 


If  no  replications  are  performed  then  Eq  (4.6)  cannot  be 
evaluated  and  Eq  (4.5)  must  use  MS^  instead  of  MSfi. 

The  calculated  P  statistics  (F^  and  F,)  are  then 

p  L 

compared  to  tabulated  P  values  to  test  the  hypothesis  of 
significant  variables  and  overall  model.  It  should  be  noted 
that  the  hypothesis  tested  is  whether  the  parameters  of  the 
variables  are  significant  from  zero.  This  is  basically 
asking  the  question  "Is  the  proposed  function  /  better  than 
no  model  at  all?”  Box  and  Draper  (1:275-280)  point  out  that 
more  than  this  simple  satistical  significance  of  the 
parameters  is  required  to  ensure  that  the  estimated  model  is 
reasonable.  It  must  be  validated  that  the  change  in  the  true 
function  over  the  design  region  is  larger  than  the  error  of 
the  estimated  model  in  the  region.  Box  and  Draper  develop  a 
conservative  criteria  to  ensure  that  this  condition  is  met. 
They  assert  that  if  the  calculated  F  is  10  times  larger  than 
the  tabulated  F  value  the  change  in  the  true  function  is 
larger  than  the  error  of  the  estimated  model.  Some 
practitioners  feel  this  is  too  conservative,  but  most  all 
tests  in  this  investigation  met  the  10  times  F  rule. 

The  final  result  of  linear  regression  is  that 
insignificant  variables  are  removed  in  a  screening  type 
process  and  that  a  parsimonious  polynomial  is  found  that 
approximates  the  investigated  relationship  over  the  restricted 
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region  of  interest.  The  next  section  explains  how  to  create 
the  experimental  design  matrix,  Z,  used  to  calculate  the 
equation  parameters. 


Fractional  Designs 

In  order  to  generate  the  y  vector  in  Eq  (4.4)  the 
mathematical  model  of  the  CSDL-I  needs  to  be  exercised  at 
specific  levels  of  the  chosen  input  variables.  The  CSDL-I 
variables  and  levels  chosen  are  described  in  the  next 
chapter.  Fractional  factorial  orthogonal  designs  were  used 
in  order  to  efficiently  generate  the  data  set.  2  q 
fractional  designs  were  implemented  on  the  input  variables, 
and  3^~q  fractional  designs  were  implemented  on  the  noise 
factors.  The  noise  factors  or  sensor  errors  are  discussed 
further  in  the  next  section  and  the  following  chapter. 

To  simplify  calculations  and  to  remove  variable  scale 
effects  the  input  variables  are  often  coded  to  values  of 
-1,0,1  for  the  low,  medium,  and  high  settings  of  the  input 
variables.  This  coding  is  accomplished  by  a  simple  linear 
transformation. 


5  5avg 
c  1/2(2*) 


(4.11) 


where 


*  =  coded  values  -1,  0,  1 
c 

*  =  uncoded  values 

^avg  “  ^^max  ^min^^ 

^  ”  ^max  ^min 
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To  illustrate  the  concept  of  fractional  designs  consider  a 
coded  design  matrix  Z  for  two  input  variables.  The  design 

2 

for  investigating  all  possible  combinations  would  require  2 
runs  and  be  of  the  following  form. 


If  the  proposed  functional  form  was 


y  =  0O  +  *  e25c2  +  612?c15c2  (4.13) 

the  design  matrix  would  need  additional  columns  to  estimate 
00  and  e12  using  Eq  (4.4).  The  complete  full  factorial  design 
would  be 

*d  ^c2  ^cl^c2 
1-1-1  1 

1  1  -1  -1  (4.14) 

1-11  -1 

ill  1 


noticing  the  mean,  0Q ,  is  estimated  by  adding  a  column  of 
ones.  The  design  matrix  in  Eq  (4.14)  is  composed  of 
orthogonal  columns.  The  dot  product  between  each  column  is 
zero.  This  orthogonal  property  allows  each  0^  to  be 
estimated  separately  without  affecting  the  values  of  the 
other  parameters.  If  $2  was  found  to  be  statistically 
insignificant  the  proposed  model  would  become 

»  *  80  ♦  ®l«cl  ♦  ®12«el«c2  <4-15> 
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The  original  values  of  6q  , 0,  and  0^2  would  remain  unchanged 
after  deleting  02*  This  is  not  true  for  nonorthogonal 
designs,  and  is  a  major  benefit  in  screening  variables. 
Another  benefit  of  orthogonal  designs  is  that  the  variance  of 
the  estimate  of  parameter  0  is  minimised.  This  is  further 
explained  in  Box  and  Draper  (1:76-80). 

If  the  proposed  function  had  three  variables  €ci'  £c2 

3 

and  £c3  a  full  factorial  design  would  require  2  or  8  runs. 

If  is  possible  to  reduce  the  number  of  runs  in  half  by 

3-1 

blocking  or  creating  a  fractional  2  design.  Again  using 
the  design  matrix  in  Eq  (4.13)  the  additional  variable  £c3 
could  be  assigned  to  the  interaction  column  T^e 

actual  experiment  would  have  at  the  low  level,  £2  at  the 
low  level,  and  £3  at  the  high  level  for  run  number  one.  The 
price  paid  for  decreasing  the  number  of  runs  from  8  to  4  is 
that  certain  interaction  terms  are  now  aliased  with  other 
effects  and  can't  be  estimated  separately.  The  most  obvious 
alias  is  that  interaction  i®  now  aliased  with  main 

effect  £c3.  The  major  assumption  is  that  higher  order 
effects  (i.e.  €ci£c2)  are  smaller  than  main  effects,  and 

therefore  affect  the  estimate  of  03  very  little.  All 
aliasing  patterns  can  be  conveniently  calculated  using  a 
generating  function.  In  the  given  example  the  function  is 
3*12.  By  applying  the  rule  that  3*3=1,  1*1=1,  etc  to  the 
generating  function  all  alias  patterns  can  be  found.  For 
example  the  following  are  present:  1*123  (mean  is  aliased 
with  interaction  1,2,3),  13*2  (interaction  13  aliased  with 


main  effect  2),  and  23#1  (interaction  23  alasied  with  main 

effect  1).  All  fractional  2*  q  designs  used  in  this 

investigation  were  standard  designs  listed  in  Box  and  Draper 

(1:164-165)  of  resolution  III  or  higher.  The  resolution  III 

simply  refers  to  the  fact  that  no  main  effects  (ie 

5cl,?c2'etc^  are  al*ase<*  wit*1  other  main  effects.  A 

resolution  III  design  is  the  minimum  required  for  estimating 

1st  order  models.  When  estimating  higher  order  models  full 

factorial  experiments  were  used  so  that  aliasing  was  not  a 

k— o 

concern.  The  fractional  3  q  design  used  in  the  noise  matrix 
was  a  resolution  III  design  taken  from  a  National  Bureau  of 
Standards  document  (3a: 13). 

The  next  section  develops  the  signal -to-noise  ratio 
which  was  used  in  Eq  (4.4)  to  investigate  the  system 
robustness  to  noise. 

Signal -to-Noise  Concepts 

One  of  the  main  objectives  of  this  investigation  was  to 
determine  and  demonstrate  methods  to  design  a  decoupled 
controller  to  be  robust  to  errors  in  specification  of  model 
parameters  and  sensor  errors.  Japanese  manufacturing 
techniques  have  proven  themselves  to  be  very  effective  in 
designing  systems  to  be  robust  to  manufacturing  variances. 
These  Japanese  statistical  tools  were  used  in  this  study  to 
help  design  the  decoupled  controller  to  be  robust  to  sensor 
errors.  The  specific  tools  used  were  the  signal-to-noise 
ratios  and  parameter  design  techniques  developed  by  Oenichi 
Taguchi  (10:167-205). 


The  first  step  in  signal -to-noise  analysis  is  the 
identification  of  design  parameters  that  are  easily 
regulated.  The  next  step  is  to  identify  factors  that  are 
difficult  or  expensive  to  regulate, i .e.  noise  factors.  In 
the  CSDL-I  example  the  fundamental  design  parameters  that 
were  easily  regulated  were  controller  gain  and  the 
arrangement  of  modes  within  the  decoupled  controller.  The 
noise  factors  identified  were  model  parameter  errors  and 
errors  in  the  sensors.  The  actual  factors  selected  could  be 
numerous  and  more  discussion  on  this  is  provided  in  the 
Choosing  Metrics  Section. 

The  next  step  is  to  construct  efficient  experiments 
(orthogonal)  that  investigate  the  effect  of  different  design 
parameters  on  the  output  for  a  range  of  noise  settings.  The 
total  experiment  consists  of  orthogonal  settings  for  the 
design  parameters  and  orthogonal  settings  for  the  noise 
factors.  For  each  design  point  associated  with  the  design 
parameters  there  is  a  corresponding  complete  set  of  orthogonal 
noise  factors.  As  a  side  point  one  should  consider  that  this 
experiment  requires  the  ability  to  enact  some  type  of  control 
over  the  noise  factors.  This  controllability  could  tempt  an 
experimenter  to  treat  the  noise  factors  as  design  parameters. 
If  the  noise  factors  were  treated  as  design  parameters  they 
would  have  to  be  maintained  at  desirable  levels  throughout 
the  life  span  of  the  system.  The  fundamental  concept  is  that 
it  is  often  easier  and  more  economical  to  set  other  design 
parameters  to  such  levels  that  the  desired  output  is  affected 


little  by  changes  in  the  noise  factors.  In  the  CSDL-I 
experiment  it  maybe  possible  to  build  sensors  accurate  to 
-  .0001  in  displacement,  but  a  sensor  accurate  only  to  -  .01 
in  maybe  cheaper,  lighter  and  more  durable.  If  the  control 
gain  and  decoupled  controller  arrangement  can  be  chosen  so 
that  the  -.01  sensor  affects  the  output  very  little  the  new 
design  will  be  cheaper,  lighter,  more  durable,  yet  with 
similar  output.  This  type  of  design  is  also  beneficial  for 
space  operation  since  the  system  degradation  would  occur  in  a 
graceful  manner.  This  is  very  important  in  the  no 
maintenance  environment  of  space. 

Once  the  experiment  has  been  conceived  and  executed,  the 
best  settings  for  the  design  parameters  can  be  obtained  by 
maximizing  a  function  derived  through  linear  regression  of 
the  observed  response.  The  key  in  applying  signal -to-noise 
techniques  is  that  the  output  variable  to  maximize  should 
reflect  not  only  robustness  to  the  noise  setting,  but  also 
reflect  system  performance.  In  the  CSDL-I  dynamic  feedback 
control  the  most  robust  system  to  senor  noise  would  be  a 
controller  with  zero  gain.  The  controller  would  not  be 
sending  any  signal  to  the  actuators  so  the  inputs  to  the 
system  would  have  no  effect.  Obviously  this  type  of  design 
is  of  little  use.  This  leads  to  the  last  step  in 
signal-to-noise  concepts,  the  signal-to-noise  ratio. 

Taguchi  (10:172)  develops  a  loss  function  to  describe  how 
deviations  from  a  target  value  contribute  to  loss  of  the 
final  design's  value.  The  loss  function  proposed  by  Taguchi 
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is  a  quadratic  function. 

l(y)  =  C(V  -  r)2  (4.16) 

where 

C  =  constant 
y  =  output 

x  =  target  value  of  y  (desired  value) 

£(«)  =  loss  function 

The  actual  loss  due  to  a  particular  design  would  be  the 
expected  value  of  Eq  (4.16). 

E U(v>)  *  C«E[(V  *  t)2]  (4.17) 

In  the  CSDL-I  the  performance  was  measured  as  a  function  of 
the  response  time  required  to  dampen  the  system.  The  target 
value  (t)  of  the  response  time  would  be  xero  since  the 
desired  performance  would  be  to  immediately  dampen  the  system 
(obviously  impractical).  In  this  case  Eq  (4.17)  reduces  to 

E(£(V))  =  CoE[(V)  ]2  (4.18) 

The  maximum  likelihood  estimate  of  the  loss  is  simply 
nn 

E U(V))  =  J  Vi/*"  (4.19) 

i  =  l 

where 

nn  =  number  of  noise  settings  (i.e.  number  of  data 
points  at  one  design  setting) 

Taguchi  then  recommends  a  transform  of  Eq  (4.19)  to  the 
following: 

nn 

SN  =  -lOlog  X 

10  4, 

i=l 
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where 


SN  =  signal-to-noise  ratio 

The  negative  sign  allows  the  loss  to  be  minimized  when  the  SN 
ratio  is  the  largest  (i.e.  less  negative);  a  convention 
adopted  in  the  literature.  The  actual  SN  ratio  used  in  this 
investigation  was: 


SN  =  -10 1 og 

to 


nn 


i  =  l 


(4.21) 


The  +101ogiQ(nn)  constant  was  dropped  for  ease  in  calculation 
since  all  noise  settings  involved  27  runs.  This  would  not 
be  done  if  the  number  of  noise  settings  changed  within  the 
experiment . 

The  next  chapter  utilizes  these  statistical  tools  to 
investigate  the  behavior  of  the  CSDL-I  model  as  developed  in 
Chapters  2  and  3.  These  tools  allow  the  decoupled 
controller's  design  parameters  to  be  set  at  values  so  as  to 
make  the  system  robust  to  input  noises. 


V.  Investigation 


Program  Development 

In  order  to  perform  the  numerical  experiments  necessary 
to  design  a  robust  decoupled  controller  a  computer  program 
was  written.  A  script  file  program  in  Pro-MATLAB  was  written 
to  enact  the  equation  of  motion  described  in  Chapter  3.  A 
documented  listing  is  provided  in  Appendix  F.  The  program 
reads  in  an  input  file  consisting  of  the  modal  transform 
matrix  U,  a  vector  V  of  the  natural  frequencies  squared,  a 
matrix  n  consisting  of  possible  decoupled  controller 
arrangements,  a  matrix  D  representing  the  global  coordinate 
input  matrix,  a  matrix  HOISE  of  sensor  noise  gains,  a  matrix 
Qc  of  the  Kalman  filter  state  weighting  values,  and  a  matrix 
Rc  of  the  Kalman  filter  control  weighting  values.  The 
outputs  are  varied,  but  the  main  ones  include  a  time  history 
of  the  state  vector  x  to  an  impulse  input,  and  the  response 
time  to  an  impulse  input  for  each  unique  setting  of  the  input 
variables.  The  time  history  of  motion  for  the  CSDL-I  was 
generated  by  applying  an  impulse  in  the  equations  of  motion, 
Eqs  (3.19)  and  (3.20).  (The  Ae  matrix  in  Eq  (3.19)  being 
specified  in  Eq  (3.65)  for  decoupled  control  and  Eq  (3.17) 
for  full -order  control.)  The  program  calculates  this  time 
history  using  a  Pro-MATLAB  subroutine,  impulse.  The  program 
was  verified  for  accuracy  by  comparing  the  eigenvalues  of  the 
close  loop  state-space  plant  matrix  Ae  to  output  for  a 
similar  program  developed  by  Robert  Calico,  PhD,  Air  Force 


Institute  of  Technology,  Wright-  Patterson  AFB,  OH.  The  new 
program  was  written  to  allow  easier  external  access  to  design 
parameters  and  more  readily  handle  statistical  analysis.  The 
following  sections  will  first  discuss  how  the  different  types 
of  metrics  were  choosen,  then  examine  the  actual  numerical 
experiment . 

Choosing  Metrics 

In  choosing  the  proper  performance  metrics,  three 
important  questions  raised  by  Box  and  Draper  (1:4)  should  be 
answered.  Which  input  variables  should  be  studied  and  in 
what  form?  How  should  response  be  measured  and  in  what  form? 
What  should  the  design  region  be?  In  trying  to  answer  these 
questions  it  becomes  apparent  that  the  possible  number  of 
variables  is  almost  unbounded.  The  possible  number  of 
variables  can  be  better  visualized  by  examining  the 
outscoping  diagram  of  the  problem.  Figure  5-1.  It  can  be 
seen  that  even  decoupled  control  becomes  a  variable  when 
examining  the  more  fundamental  question  of  maintaining  the 
alignment  of  a  space  structure.  The  problem  was  immediately 
narrowed  down  to  the  investigating  of  how  modal  arrangements 
within  the  decoupled  controller  and  controller  gain  affect 
the  system  robustness.  Once  the  problem  was  narrowed  down  to 
a  size  manageable  for  allotted  resources  specific  inputs 
could  be  chosen.  The  following  subsections  will  discuss 
noise  inputs,  then  response  measures  and  finally  design 
parameters. 


Figure  5-1  Outsooping  of  Problem 


Input  Noise  Variables.  Two  types  of  noise  inputs  were 
selected.  The  first  type  of  noise  was  modelling  error.  Poor 
specification  of  the  model  could  take  on  many  forms.  A 
structural  model  often  contains  error  due  to  the 
approximation  of  a  continuous  structure  as  a  finite  element 
model.  This  type  of  error  results  in  higher  natural 
frequencies  being  wrong  and  some  completely  missing.  This 
type  of  noise  was  simulated  by  only  designing  the  controllers 
based  on  the  first  eight  natural  frequencies.  However,  time 
performance  histories  were  based  upon  the  total  model  of 
twelve  natural  frequencies.  A  robust  design  should  be  able 
to  enact  good  control  on  a  structure  even  though  the  design 


5-3 


was  based  on  incomplete  knowledge.  This  first  type  of  noise 
was  always  present  and  was  not  varied  in  a  signal -to-noise 
type  experiment,  but  was  treated  as  natural  or  pure 
experimental  error. 

The  second  type  of  noise  input  selected  was  incorrect 

sensor  readings.  A  sensor  may  not  be  calibrated  correctly 

and  give  incorrect  information  to  the  feedback  system.  A 

robust  design  should  still  be  able  to  enact  adequate  control 

on  a  structure  even  with  slightly  wrong  information.  Since 

only  position  sensors  were  used,  the  errors  could  only  be 

induced  in  position  readings.  Velocity  or  acceleration 

sensors  would  have  different  noise  effects  than  position 

errors.  Th^  position  errors  should  be  a  little  more  serious 

than  velocity  errors  since  the  estimator  differentiates  the 

position  readings  to  derive  the  other  states.  The  errors 

were  induced  by  multiplying  the  actual  sensor  readings  by 

coefficients  of  .8,  1.0  and  1.2.  The  .8  coefficient  causes 

the  position  sensor  to  always  read  low  by  20%.  It  is 

important  to  note  that  the  errors  were  supplied  in  the  local 

coordinate  frame  i  and  then  transformed  into  the  appropriate 

perturbed  C  matrix  using  the  rotation  matrices  in  Eq  (2.13) 

(2.14)  and  (2.15).  The  Kalman  estimator  gain  design  was 

based  upon  noise  being  decoupled  in  the  state-space 

coordinate  frame;  a  common  assumption  since  the  actual  noise 

distribution  is  usually  unknown.  A  robust  design  should  be 

insensitive  to  real  noises.  This  second  type  of  noise  was 

6-3 

varied  in  a  signal -to-noise  type  experiment  using  a  3 
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fractional  design. 


Response  Measures.  The  final  desired  output  was 
previously  defined  in  Chapter  2  as  maintaining  the  line  of 
sight  (LOS).  The  question  is  how  to  measure  the  LOS.  This 
is  a  difficult  question  to  answer.  The  CSDL-I  model  was 
small  enough  that  it  allowed  performing  various  experiments 
on  different  output  variables  for  various  settings  of  input 
variables.  The  criterion  for  a  good  metric  for  LOS  are; 
ability  to  measure  the  deviation  of  node  one  from  it’s 
undisplaced  position,  and  ability  to  measure  the  time 
required  to  reduce  deviation.  Also  the  metric  itself  must 
have  a  large  value  spread  for  differing  inputs. 

The  deviation  of  node  one  was  measured  by  calculating 
the  length  and  angle  of  a  position  vector  from  the 
undisplaced  location  for  the  time  history  of  node  one.  A 
representative  trace  of  the  displacement  is  shown  in  Figure 
5-2.  The  following  derived  metrics  were  investigated: 
average  position  vector  length  (radius)  over  a  20  second 
period,  average  angle  over  a  20  second  period,  standard 
deviation  of  radius,  standard  deviation  of  angle,  a  response 
time  to  a  minimum  moving  average,  and  an  average  velocity  of 
node  one.  The  final  performance  metric  based  on  the  above 
mentioned  criteria  was  the  r  onse  time.  The  response  time 
(tr)  is  the  time  when  the  average  radial  displacement  (using 
a  1.48  sec  window)  remained  below  .05  displacement  units.  A 
plot  of  210  different  subcontroller  arrangements  is  shown  in 
Figure  5-3.  The  plot  shows  a  good  spread  of  data  from  6  to 
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20  seconds.  The  tr  metric  is  a  direct  measure  of  the  amount 
of  time  to  reduce  node  one  displacement,  and  also  provides  a 
measure  of  the  magnitude  of  node  one  displacement  since  it  is 
a  computed  time  based  on  an  average  radial  displacement. 


Figure  5-2  LOS  Trace  for  Arrangement  123-45-678 


The  response  time  was  determined  as  an  appropriate 
output  variable,  but  in  order  to  design  a  robust  controller  a 
signal -to-noise  (SN)  ratio  variable  was  necessary.  By 
applying  Eq  (4.20)  to  tr,  a  signal -to-noise  ratio  was 
calculated.  All  regressions  in  this  investigation  used  both 
tr  and  SN  of  tr  as  output  responses. 

Design  Parameters.  In  choosing  design  parameters  it  was 
decided  to  only  investigate  the  gain  associated  with  the 
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Figure  5-3  Response  Time  vs  Rank 


controller,  Kc  in  Eq  (3.3),  and  not  the  estimator,  Kf  in  Eq 

(3.4).  In  exploring  the  effects  of  controller  and  estimator 
2 

gains  a  small  3  design  was  performed  on  a  few  specific 
subcontroller  arrangements.  The  input  variables  were  Qq 
matrix  for  the  estimator,  and  the  Qq  matrix  for  the 
controller,  Eq  (3.8).  As  stated  before  the  Qq  matrix  was 
always  a  diagonal  matrix.  In  this  initial  exploratory 
screening  all  the  diagonal  elements  were  changed 
simultaneously  for  the  estimator  denoting  this  as  a  single 
variable  Q0^»  and  for  the  controller  denoted  as  Qoc* 

A  regression  on  the  SN  ratio  was  performed  and  the  results  are 
in  Tables  0-1  and  0-2,  Appendix  G.  A  graph  of  the  contours 
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associated  with  constant  values  of  the  tr  SN  ratio  are  shown 
in  Figure  5-4.  The  results,  show  as  expected,  that  the 
larger  total  compensator  gain  gives  the  largest  SN  ratio. 


Figure  5-4  Estimator  vs  Controller  State  Weighting 

It  was  originally  hoped  that  Figure  5-4  would  show  a  clear 
stationary  ridge.  If  this  had  been  the  case,  it  could  have 
been  concluded  that  the  estimator  and  controller  gain  were 
only  measures  of  a  more  fundamental  variable,  namely 
compensator  gain.  However,  this  was  not  evident,  but  in 
order  to  reduce  the  number  of  inputs,  was  fixed  at  a 
coded  1  (uncoded  value  of  20),  to  coincide  with  the  maximum 
signal -to-noise  ratio.  In  all  subsequent  experiments  was 
20.  Also  RQ£  was  set  to  an  identity  matrix  for  reasons 
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explained  below. 


After  it  was  decided  to  investigate  the  controller  gain 

effects,  Kc,  the  inputs  to  Kc  were  further  restricted  by  only 

varying  Qq,  and  not  Rq.  A  small  3  experiment  was  performed 

varying  the  diagonal  elements  of  Qq  and  Rq  denoting  these  as 

single  variables  Q _  and  R^  .  The  regression  analysis  on  the 

oc  oc 

tr  SN  ratio  is  provided  in  Tables  G-2  and  G-3,  Appendix  G. 
Again  the  contours  for  constant  values  of  SN  ratio  are 
plotted  in  Pigure  5-5. 


Figure  5-5  Control  vs  State  Weighting 


It  was  expected  that  the  contours  would  show  a  stationary 
ridge  since  the  controller  gain  Kc  is  directly  a  function  of 
R  and  Q  as  described  by  Egs  (3.9)  and  (3.10).  Again  this 

OC  OC 
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was  not  evident,  but  in  order  to  reduce  the  number  of  inputs, 

R  was  restricted  to  an  identity  matrix.  While  the 
oc 

regression  analysis  indicated  that  the  effects  of  Q  .  and  R _ 

oc  oc 

are  very  similar,  Q  was  chosen  as  a  design  parameter  because 

o 

its  matrix  dimension  is  larger  than  the  RQ  matrix.  The 
larger  dimension  allows  for  a  more  detailed  adjustment  on  the 
states  associated  with  the  natural  frequencies. 

It  is  interesting  to  note  that  by  examining  Eq  (3.8)  the 
above  results  could  have  been  deduced  for  the  response  time, 
but  not  necessarily  the  tr  SN  ratio  (the  response  time  is  a 
function  of  the  size  of  gain).  In  Eq  (3.8)  it  can  be  seen 
that  a  small  Rq  compared  to  a  large  Qq  simply  means  that  the 
most  performance  improvement  comes  from  reducing  the  amount 
of  time  to  dampen  the  structure.  The  size  of  the  control 
effort  becomes  a  less  important  factor.  This  is  the  the  same 
result  (i.e.  small  Rq,  large  Qq)  that  was  experimentally 
determined  for  the  tr  SN  ratio.  While  this  investigation 
didn't  treat  actuator  energy  as  a  response  metric  it  is 
important  and  should  be  mentioned.  If  actuator  energy  had 
<  _en  chosen  as  an  output,  a  strictly  small  RQ  and  large  Qq 
would  not  have  been  the  best  setting.  A  good  range  of  values 
for  Rq  and  Qq  could  have  been  selected  by  using  both  tr  and 
actuator  energy  in  an  initial  screening  type  experiment.  If 
a  second  response  metric  of  total  actuator  or  maximum 
actuator  energy  contour  plot  was  overlaid  on  Figure  5-5,  a 
region  defining  appropriate  values  of  RQ  and  Qq  would  have 
been  defined.  This  region  would  represent  the  compromise 


between  quickly  dampening  system  and  the  actuator  energy. 

This  side-point  emphasizes  the  fact  that  simply  maximizing 
the  relationship  shown  in  Pigure  5-5  would  not  always  be  a 
wise  decision  since  extremely  large  Qq  would  over  burden  the 
actuators . 

In  summary,  the  response  time,  tr,  and  signal -to-noise 
ratio  of  tr  were  chosen  as  the  response  metrics.  The  noise 
inputs  were  sensor  noise  variation  and  modal  truncation  of 
the  highest  four  modes.  The  design  parameters  included  the 
16  diagonal  components  of  the  Qq  matrix  used  in  calculating 
the  controller  gain  matrix  Kc,  and  the  particular  arrangement 
of  the  lowest  eight  modes  into  a  3-2-3  decoupled  controller. 
While  this  seems  like  a  narrow  slice  of  the  total  design 
problem,  the  possible  combinations  required  to  investigate  is 
still  quite  astounding.  This  choice  of  variables  represents 
a  possible  3®  x  21®  x  560  or  2.7(10®)  combinations.  Add  to 
this  the  fact,  the  mathematical  modeling  requires  three 
possible  representative  impulses,  the  design  space  grows  to 
8.03(10®)  combinations.  The  CSDL-I  program  written  for  this 
investigation  required  approximately  30  seconds  CPU  time  to 
generate  the  time  history  necessary  to  compute  tr  and  tr  SN 
ratio.  This  equates  to  76,354  years  of  CPU  time  to  generate 
a  complete  data  set.  The  power  of  RSM  is  utilized  in  the 
next  section  to  intelligently  probe  and  screen  the  design 
space  to  arrive  at  a  good  design. 

Experiment 


In  order  to  find  the  best  combinations  of  design 


parameters  out  of  the  enormous  possibilities,  two  functions 
must  be  accomplished;  screen  design  parameters,  and  develop  a 
relationship  for  only  the  important  parameters.  The 
following  subsections  will  first  examine  the  subcontroller 
mode  arrangement  then  the  controller  gain. 

Mode  Arrangements.  The  first  parameter  screened  was  the 
arrangement  of  the  eight  highest  modes  in  the  3-2-3 
controller.  Traditional  RSM  screening  requires  continuous 
independent  variables.  This  makes  it  impossible  to 
investigate  each  subcontroller  separately  as  an  independent 
parameter.  For  example,  when  different  modes  are  assigned  to 
subcontroller  one,  another  subcontroller's  modes  must  also 
change.  This  means  the  assignment  of  modes  to  subcontrollers 
must  be  treated  as  a  single  qualitative  parameter  consisting 
of  560  possible  combinations.  While  it  would  be  possible  to 
investigate  all  560  combinations  for  the  CSDL-I,  a  larger 
more  realistic  model  would  require  too  much  computational 
effort.  This  lead  to  the  development  of  a  heuristic  to 
screen  the  design  space  without  investigating  all  560  design 
points . 

In  decoupled  control  the  individual  subcontrollers  don't 
exercise  the  same  amount  of  control  over  their  particular 
assigned  modes.  Some  of  the  subcontrollers  must  expend 
actuator  energy  to  maintain  orthogonality  between 
subcontrol lers  so  as  to  not  destabilize  each  other.  The 
condition  used  in  Eq  (3.34)  enforces  maximum  controllability 
in  subcontroller  one,  and  minimum  controllability  in 
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subcontroller  three.  Likewise,  it  enforces  maximum 
observability  in  subcontroller  three  and  minimum 
observability  in  subcontroller  one.  It  is  also  known  that 
the  lower  frequency  modes  are  more  effective  in  controlling 
the  structure  than  the  higher  modes.  This  logic  leads  to 
placing  the  lowest  modes  (i.e.  1,2,3)  into  subcontroller  one, 
the  next  lowest  (i.e.  4,5)  into  subcontroller  two,  etc.  This 
ascending  order  (AO)  design  is  usually  a  good  design  based  on 
traditional  practices.  The  important  information  it  doesn't 
use  is  that  certain  modes  tend  to  be  aligned  and  prefer  to  be 
controlled  together  rather  than  in  separated  subcontrollers. 
Calico  (3:43).  The  following  rules  determine  a  controller 
arrangement  that  keeps  modes  that  are  aligned  together 
grouped  in  the  same  subcontroller. 

Arrangement  Heuristic.  These  rules  are  also  shown 
in  a  flow  chart  type  format  in  Figure  5-6. 

1)  Always  place  mode  one  in  subcontroller  one. 

note:  This  rule  is  based  on  the  fact  that  it  is 
desirable  to  place  the  most  controllable  mode  into  the 
subcontroller  able  to  enforce  the  most  control. 

2)  Place  the  lowest  modes  into  subcontrollers  one  (i.e. 
1,2,3).  Generate  all  possible  combinations  in 
subcontroller  two  while  keeping  subcontroller  one  modes 
fixed.  Exercise  model  for  a  single  impulse  input  and 
select  the  arrangements  with  the  two  shortest  time  of 
responses.  For  book  keeping  purposes  designate  the  two 
controller  arrangements  as  2A  and  2B.  The  2A 
arrangement  having  the  shorter  response  time. 

note:  See  note  for  rule  3,  page  5-14. 
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3)  Place  the  highest  modes  into  subcontroller  one  while 
still  observing  rule  one  (i.e.  1,7,8).  Generate  all 
possible  combinations  in  subcontroller  two  while  keeping 
subcontroller  one  modes  fixed.  Exercise  model  for  the 
same  impulse  input  as  rule  2  and  select  the  arrangements 
with  the  two  shortest  time  of  responses.  Por  book 
keeping  purposes  designate  the  two  controller 
arrangements  as  3A  and  3B.  The  3A  arrangement  having 
the  shorter  response  time. 

note:  The  above  two  rules  investigate  the  design  space 
at  two  extremes.  In  rule  two  the  most  controllable 
modes  are  placed  in  subcontroller  one  and  in  three  two 
the  least  controllable  modes  are  assigned  to 
subcontroller  one. 

4)  If  2A  and  3A  contain  different  mode  groups  in 
subcontroller  two  fix  these  subcontroller  two 
combinations  and  generate  all  possible  arrangements  for 
subcontroller  one  and  three.  If  2A  and  3A  contain  the 
same  mode  groups  in  subcontroller  two  replace  the  longer 
response  time  arrangement,  either  2A  or  3A,  with  2B  or 
3B  respectively.  Generate  all  possible  combinations  in 
subcontroller  one  and  three.  Exercise  the  model  for  a 
single  impulse  input  and  select  the  three  shortest  time 
responses  for  each  unique  controller  two  combination. 

note:  The  above  rule  has  made  use  of  the  fact  that  the 
two  modes  in  subcontroller  two  had  superior  performance 
at  the  extreme  ends  of  the  design  space.  The  modes 
prefer  to  be  in  the  same  controller. 


5)  Take  the  six  arrangements  selected  from  rule  4  and 
exercise  the  model  for  a  different  impulse  direction. 
Select  the  three  best  arrangements  from  this  second 
impulse . 


note:  The  above  rule  tries  to  select  arrangements  that 
are  robust  to  impulse  direction.  While  impulse 
direction  wasn't  a  noise  variable  in  this  investigation 
rule  5  is  essential  to  selecting  a  good  arrangement.  It 


is  possible  that  an  arrangement  may  have  superior 
performance  to  a  particular  impulse  by  supplying  control 
energy  to  a  single  mode  that  is  especially  sensitive  to 
that  particular  impulse.  An  impulse  from  a  different 
direction  will  reveal  this  situation  as  the  controller 
arrangement  will  have  significantly  degraded  performance 
for  the  different  impulse. 

The  above  heuristic  was  applied  to  the  CSDL-I  model  yielding 
top  three  arrangements  of  125-47-368,  126-47-358,  and 
158-47-236.  A  flow  chart  showing  this  process  for  the  CSDL-I 
example  is  given  in  Figure  5-7.  The  heuristic  required 
exercising  the  model  a  total  of  43  times.  A  full  factorial 
experiment  was  also  performed  on  the  model  verifing  the 
heuristic . 

It  would  seem  logical  to  select  only  one  best 
arrangement  from  the  heuristic,  but  this  could  lead  to  some 
problems.  The  CSDL-I  model  was  exercised  with  all  parameters 
fixed  except  the  subcontroller  arrangement.  It  could  be 
possible  that  there  are  interaction  effects  between 
controller  gain  and  subcontroller  arrangement.  By  carrying 
forward  at  least  three  possible  arrangements  into  the  next 
screening  experiment  these  second  order  effects  may  be 
revealed,  but  there  is  no  guarantee.  The  final  design  cannot 
be  called  the  "best"  possible  design,  but  only  a  "good" 
design . 


goto  goto  goto 
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Pigure  5-6  Modal  Arrangement  Heuristic 
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Figure  5-6  Modal  Arrangement  Heuristic  (continued) 


Figure  5-7  Modal  Arrangement  Heuristic  Example 
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Modal  Arrangement  Heuristic  Example  (continued) 


Controller  Gains.  The  next  design  parameters  screened 


were  the  16  diagonal  components  of  the  matrix  used  to 
calculate  controller  gains,  Kc,  as  defined  in  Eq  (3.8). 

A  group  screening  experiment  was  performed  on  the  16 

factors.  A  group  screening  experiment  logically  arranges  the 

individual  factors  into  groups  and  then  varies  the  groups  as 

a  whole  to  determine  which  groups  are  statistically 

significant.  A  basic  assumption  is  that  the  important 

individual  factors  will  be  contained  in  the  significant 

groups.  The  logical  grouping  of  variables  consisted  of 

sorting  the  Qq  elements  by  subcontroller  and  by  velocity  or 

position  state.  The  six  groups  were  varied  from  20  for  a  low 

6- 2 

value  to  40  for  a  high  value  using  a  2jy  fractional 
experiment  The  experimental  design  matrix  taken  from  Box  and 
Draper  (1:164)  is  listed  in  Appendix  G.  Three  output 
variables  were  calculated,  time  of  response,  variance  of  time 
of  response,  and  signal-to-noise  ratio  for  time  of  response. 

6-3 

The  standard  3  fractional  design  previously  mentioned  was 
used  for  the  noise  variable.  The  experiment  was  run  for  all 
three  possible  arrangements  and  the  summary  of  results  is 
contained  in  Tables  G-5  and  G-6,  Appendix  G  for  arrangement 
126-47-358.  This  experiment  was  used  to  identify  important 
groupings  of  the  16  Qo  elements  and  to  identify  which  of  the 
three  arrangements  to  carry  forward  for  further 
experimentation . 

Arrangement  126-47-358  was  chosen  as  the  best  overall 
arrangement  based  on  the  average  SN  ratio  and  the  average 
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variance  of  the  time  response.  Arrangement  125-47-368  had  an 
average  SN  ratio  of  -28.3  and  an  average  variance  of  time 
response  of  1.043.  Arrangement  126-47-358  had  an  average  SN 
ratio  of  -29.2  and  an  average  variance  of  time  response  of 
.0467.  Arrangement  158-47-236  had  an  average  SN  ratio  of 
-29.4  and  an  average  variance  of  time  response  of  .096.  The 
actual  performance  of  the  three  arrangements  was  really  quite 
similar.  The  SN  ratio  which  was  the  primary  output  metric 
for  designing  a  robust  controller  had  a  spread  of  only  4%. 

If  the  SN  ratio  had  been  strictly  used,  arrangement  125-47-368 
would  have  been  chosen  since  it  had  the  largest  SN  ratio.  By 
examining  the  results  a  little  closer  a  better  choice  was 
made.  The  actual  time  response  of  all  three  arrangements  was 
already  determined  as  adequate  by  meeting  an  objective 
standard.  Arrangement  126-47-358  had  the  lowest  variance  of 
all  three,  therefore,  it  was  selected  as  the  best  arrangement 
based  on  its  robustness  to  sensor  error.  Its  signal-to-noise 
ratio  wasn't  superior  only  because  its  mean  response  time  was 
slightly  higher  than  the  mean,  for  arrangement  125-47-368,  but 
as  previously  mentioned,  this  was  already  adequate  and  could 
be  reduced  by  increasing  feedback  gain. 

By  examining  the  analysis  of  variance  Table  G-6  for  tr 
and  SN  in  Appendix  G  it  can  be  seen  that  the  statistically 
significant  groups  for  arrangement  126-47-358  were  the  Qq 
groups  for  velocity  state  for  modes  in  subcontroller  one  and 
the  groups  for  velocity  state  for  modes  in  subcontroller 
three.  It  seems  reasonable  that  the  velocity  states  should 
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be  the  most  significant  since  in  Eq  (2.25)  modal  damping  is 
associated  with  velocity  states.  These  two  significant 
groups  represent  six  individual  elements  of  the  Qq  matrix. 

The  six  elements  or  factors  were  then  investigated  to 

determine  which  of  the  six  were  significant,  since  the 

previous  experiment  only  determined  that  the  groupings  werj 

significant.  The  six  factors  were  varied  from  a  20  for  a  low 
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value  to  a  40  for  a  high  value  using  a  2  fractional 
experiment.  The  factors  from  the  insignificant  groups  were 
set  to  either  the  high  or  low  value,  whichever  maximized  the 
SN  ratio.  For  example,  the  three  Qq  elements  for  position 
states  of  the  modes  in  subcontroller  one  were  all  set  to  the 
low  value  of  20  since  the  SN  equation  coefficient  of  this 
group  was  -.0113.  It  is  important  to  set  all  the  individual 
factors  within  the  group  to  the  same  value  so  that  the 
effects  of  the  insignificant  groups  are  aliased  only  with  the 
mean.  Design  points  at  the  center  were  used  to  check  for  any 
quadratic  terms. 

By  examining  the  analysis  of  variance.  Table  G-7  and  G-8 
Appendix  G,  for  tr  and  SN  it  can  be  seen  that  the  only 
significant  factors  were  the  Qq  elements  for  the  velocity 
state  of  mode  3  in  subcontroller  three,  and  the  velocity 
state  of  mode  8  in  subcontroller  three.  These  were  the  only 
statistically  significant  factors,  but  velocity  state  of  mode 
2  in  subcontroller  one  was  almost  an  order  of  magnitude  more 
significant  than  any  other  factor,  therefore  it  was  carried 
forward  into  another  experiment.  It  is  interesting  to  note 


that  this  result  is  not  what  would  have  been  expected. 
Subcontroller  one  has  the  most  control  authority  available, 
thus  it  would  be  expected  that  the  modes  in  subcontroller  one 
would  be  the  significant  modes  in  regulating  the  system. 
Evidently,  the  synergistic  effect  of  modal  grouping  and  the 
particular  LOS  requirement  made  the  Qq  elements  associated 
with  the  velocity  states  of  modes  3  and  8  the  most  important. 
This  result  would  have  been  difficult  to  determine  without 
using  RSM . 

The  final  three  significant  factors,  velocity  state  mode 
3,  velocity  state  mode  8,  and  velocity  state  mode  2  were 

3 

further  investigated  for  interaction  terms  using  a  2 
full -factorial  experiment.  It  was  found  that  only  the 
velocity  state  mode  3,  velocity  state  mode  8,  and  their 
interaction  terms  were  significant  as  shown  in  Table  G-9  and 
G-10,  Appendix  G.  A  final  functional  form  was  determined  as: 

SN  =  -29.2  +  .086Qv3  +  .071Qy8  +  ,059Qv3Qv8  (5.1) 

tr  =  5.57  -  .057Qv3  -  .047Qv8  -.04Qy3Qy8  (5.2) 

where 

Qy3  =  diagonal  element  of  Qq  associated  with  velocity 
state  of  mode  3  in  subcontroller  three 
Qy8  =  diagonal  element  of  Qq  associated  with  velocity 
state  of  mode  8  in  subcontroller  three 

In  Eq  (5.1)  there  is  no  replication  of  the  design  points 
since  the  SN  ratio  is  calculated  across  all  of  the  noise 
settings.  The  absence  of  replications  prevents  the 
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computation  of  MS  (pure  error,  Eq  (4.8))  for  the  SN  ratio 

e 

regression.  In  Eq  (S.2)  there  is  replication  of  the  design 

points  associated  with  noise.  The  estimation  of  the  standard 

1/2 

deviation  is  (MSe)  which  is  .1547.  Typically,  this  large 
1/2 

value  of  (MSe)  '  would  suggest  that  the  coefficients  in  Eq 

(5.1)  are  of  little  use  because  of  such  a  large  underlying 

1/2 

error.  In  this  particular  case  the  error  that  the  (MSe) 
is  measuring  is  the  -  20%  sensor  readings  induced  by  the 
noise  matrix.  This  large  deviation  may  or  may  not  be 
representative  of  real  errors  in  a  physical  system,  thus  the 
coefficients  may  actually  be  larger  than  physical  underlying 
error.  The  important  observation  to  be  made  from  Eq  (5.1)  is 
that  the  SN  ratio  is  maximized  by  setting  Qv3  and  Qv8  to 
their  high  levels.  It  is  also  instructive  to  observe  that 
even  at  the  high  levels  of  Qy3  and  Qvg  the  final  SN  ratio  is 
changed  very  little.  The  Qq  factors  had  minimal  effect  on 
robustness  of  the  system  over  the  ranges  of  Qq  selected.  A 
different  range  for  Qq  might  yield  different  results.  Most 
of  the  increase  in  SN  ratio  was  due  to  selecting  a  robust 
arrangement . 

The  error  distribution  of  noise  in  the  sensors  can  be 
further  investigated  by  plotting  the  residuals  of  Eq  (5.2). 

A  residual  plot  is  the  graph  of  differences  between  the 
actual  data  points  and  the  predicted  values  versus  the 
predicted  values.  If  the  errors  have  constant  variance  and 
the  model  is  specified  correctly,  the  plot  should  have  a 


random  scattering  of  points.  The  residual  plot  is  shown  in 
Figure  5-8. 


It  can  be  seen  that  the  range  of  errors  becomes  smaller  as 
the  magnitude  of  the  response  increases.  This  would  suggest 
a  transformation  of  the  output  variable.  A  possible 
transform  might  be  log(tr(max)  -  tr).  This  was  not  actually 
done  since  the  range  of  differences  between  actual  data 
points  and  predicted  points  was  so  small,  but  other 
situations  may  require  such  transforms. 

At  the  end  of  the  experiment  the  final  design  determined 
consisted  of: 
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1)  A  3-2-3  controller  arrangement  with  modal  arrangement  of 
126-47-358. 

2)  An  identity  matrix  for  each  of  the  R  matrices  used  in 

o 

calculating  the  control  gains  Kc  and  estimator  gains  Kf 

3)  A  diagonal  matrix  of  20  for  each  of  the  Qq  matrices  used 
in  calculating  the  estimator  gains  Kf. 

4)  The  following  Q  matrices  used  in  calculating  the 

o 

controller  gains  Kc. 

subcontroller  one:  20  0  0  0  0  0 

0  20  0  0  0  0 

0  0  20  0  0  0 

0  0  0  20  0  0 

0  0  0  0  40  0 

0  0  0  0  0  40 


subcontroller  two:  [40  0  0  0 

0  40  0  0 

0  0  40  0 

0  0  0  40 


subcontroller  three:  40  0  0  0  0  0 

0  40  0  0  0  0 

0  0  40  0  0  0 

0  0  0  40  0  0 

0  0  0  0  40  0 

0  0  0  0  0  40 


After  the  final  design  was  determined,  a  verification  of 
the  adequacy  of  the  performance  was  accomplished.  In  the 
next  chapter  the  performance  of  this  design  is  compared  to 
other  designs. 


VI.  Results 


The  final  robust  design  as  described  in  the  previous 
chapter  was  evaluated  against  three  additional  structural 
control  conditions.  The  first  additional  condition  was  the 
structure  with  no-control  representing  the  worst  possible 
stable  decoupled  controller  performance.  The  next  condition 
was  a  full -order  controller  designed  using  total  knowledge  of 
the  system.  A  control  system  based  upon  all  twelve  modes 
using  LQR/LQE  theory  was  implemented.  All  of  the  weighting 
matrices  except  for  Qq/  used  to  compute  Kc,  were  the  same  for 
the  full -order  controller  as  for  the  decoupled  controller. 

The  Qo  matrix  was  a  diagonal  matrix  of  35  representing  the 
center  of  the  design  space  used  in  investigating  the 
decoupled  controller.  This  full -order  design  represents  the 
best  possible  decoupled  controller  performance.  The 
full -order  controller  doesn't  have  any  residual  modes  and 
doesn't  have  to  expend  actuator  energy  to  maintain 
subcontroller  orthogonality.  The  last  additional  condition 
was  an  ascending  order  decoupled  controller  based  on  the 
first  eight  modes.  The  ascending  order  design  represents  the 
traditional  design  methodology.  The  ascending  order  design 
was  constructed  using  LQR/LQE  theory  as  previously  developed 
with  the  same  weighting  matrices  except  Qq,  used  to  calculate 
Kc.  The  Qq  matrix  was  again  a  diagonal  matrix  of  35  to 
represent  the  center  of  the  design. 

In  Figure  6-1  the  moving  average  time  history  of  all 
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Figure  6-1  Moving  Average  of  Radius  vs  Time 


four  conditions  are  plotted.  It  is  clear  to  see  that  the 
uncontrol  led  CSDL-I  is  extremely  lightly  damped  as  the  moving 
average  decreases  very  little  over  the  20  second  interval. 

The  full-orcler  controller  displays  superior  performance  as 
would  be  expected.  This  model  has  perfect  knowledge  of  the 
system,  therefore,  it  is  capable  of  enacting  efficient 
actuator  authority.  The  ascending  order  design  and  robust 
design,  overall  displayed  fair  performance.  The  robust 
design  is  significantly  better  than  the  ascending  order 
design  at  long  durations.  While  one  cannot  call  the  robust 
design  the  optimum  design  due  to  numerous  assumptions 
preventing  rigorous  proof  of  optimal lity,  the  robust  design 
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can  be  considered  a  ''good"  design  based  upon  the  results  o£ 
Figure  6-1. 

In  Figure  6-2  tr  is  plotted  against  the  27  different 
sensor  noise  settings  for  the  robust  design,  ascending  order 
design,  and  the  full -order  design. 


Figure  6-2  Reponse  Time  vs  Noise  Settings 


The  no-control  design  was  obviously  not  plotted  since  it  used 
no  sensors  or  actuators.  This  figure  demonstrates  the 
robustness  of  each  design.  The  full -order  design  is 
extremely  robust  to  sensor  error  as  shown  in  the  graph.  The 
ability  to  perfectly  know  the  model  allowed  the  system  to 
drive  the  error  states  to  zero  very  quickly.  The  sensor 
error  gain  multiplied  by  a  zero  displacement  is  still  zero. 


The  robust  decoupled  design  is  also  extremely  robust  to 
sensor  noise.  In  fact,  the  only  major  disadvantage  of  the 
robust,  decoupled  design  over  a  full -order  design  based  upon 
Figure  6-1  and  6-2  is  that  the  decoupled  controller  is  slower 
in  dampening  the  structure.  The  ascending  order  design  is 
not  a  robust  design.  The  particular  settings  of  noise  caused 
wild  fluctuations  in  performance.  This  figure  shows  that 
modal  grouping  in  decoupled  control  has  a  dramatic  affect  on 
sensitivity  to  sensor  noise.  The  design  techniques  in  this 
investigation  produced  a  decoupled  controller  exhibiting  more 
robustness  than  those  of  traditional  design  practices. 

Before  deciding  to  implement  any  controller  design  the 
performance  required  of  the  physical  actuators  must  be  known. 
In  Figure  6-3  a  plot  of  the  time  history  of  leg  9  actuator 
force  is  plotted.  The  plot  was  created  from  the  state  vector 
data  for  an  impulse  input  with  no  sensor  noise  applied.  The 
decoupled  controller  design's  force  time  history  was  computed 
by  summing  the  individual  subcontroller  actuator  inputs  as 
described  in  Eq  (3.40).  The  full-order  force  time  history 
was  computed  using  Eq  (3.6).  As  expected  the  full -order 
model  is  very  efficient  in  using  actuator  energy  to  control 
the  structure  as  it  doesn't  have  to  expend  energy  to  maintain 
subcontroller  orthogonality.  It  was  surprising  to  see  that 
the  robust  design  was  inputting  significantly  more  energy 
into  the  actuators  than  the  ascending  order  design.  The  gain 
matrices  were  computed  using  essentially  the  same  weighting 
matrices.  The  different  arrangements  observe  the  system  and 
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Figure  6-3  Leg  9  Actuator  Force  Time  History 

try  to  enforce  control  on  the  system  differently,  thus  the 
LQR/LQE  theory  computes  different  Kc  and  Kf  matrices.  It  is 
probable  that  the  ascending  order  design  has  subcontroller 
actuator  inputs  that  are  canceling  each  other  out.  This 
would  also  explain  the  inability  of  the  ascending  design  to 
dampen  the  structure.  From  Figure  6-3  it  can  be  seen  that 
the  price  paid  for  decoupled  control  in  the  robust  design  is 
a  larger  actuator  force  requirement.  As  noted  earlier  if 
actuator  energy  had  been  a  response  parameter  it  may  be 
possible  to  find  a  subcontroller  arrangement  that  was  less 
robust,  but  also  used  less  actuator  authority. 


VII.  Conclusions  and  Recommendations 


In  conclusion,  the  statistical  design  methods  developed 
in  this  investigation  have  produced  a  new  decoupled 
controller  design  superior  to  a  traditional  design,  namely 
the  ascending  order  design.  The  statistical  design  mehtods 
employed  signal -to-noise  ratios,  orthogonal  numerical 
experiments,  linear  regression  and  a  modal  arrangement 
heuristic  to  efficiently  explore  the  large  design  space. 

These  methods  reduced  the  computational  efforts  from  a 
possible  76,354  years  of  CPU  time  for  a  global  search  method 
to  only  16.5  Hrs  of  CPU  time  for  the  CSDL-I  model.  All  of 
the  design  techniques  used  were  general  enough  to  apply  to 
any  size  of  structural  control.  It  should  be  noted  that  the 
modal  arrangement  heuristic  was  only  demonstrated  on  an  eight 
mode  3-2-3  subcontroller  arrangement  and  needs  to  be  further 
investigated  on  other  models,  but  should  still  be  applicable 
to  these  other  cases. 

The  first  recommendation  for  further  study  is  obviously 
to  try  the  modal  arrangement  heuristic  on  other  subcontroller 
arrangements.  The  3-2-3  arrangement  was  very  convenient 
since  the  middle  controller  contained  only  two  modes.  This 
allowed  the  screening  techniques  to  directly  generate  pairs 
of  modes  that  preferred  to  be  together.  For  a  larger 
subcontroller  system,  say  20-15-20,  the  exact  procedure 
should  yield  good  results,  except  the  screening  techniques 
would  generate  a  set  of  15  modes  that  prefer  to  be  grouped 


together.  For  a  larger  and  more  divided  subconlroller 
system,  say  30-15-10-30,  the  heuristic  should  still  apply  by 
simply  grouping  the  30  and  15  mode  subcontrollers  together 
and  treating  them  as  the  first  subcontroller.  The  10  and  30 
mode  subcontroller  would  then  be  treated  as  the  second  and 
third  subcontroller,  respectively.  Since  the  heuristic  isn't 
based  upon  strict  application  of  theory  these  possibilities 
should  be  exercised  on  other  models  and  compared  to  the  full 
factorial  results. 

The  second  recommendation  for  further  study  is  to 

develop  a  response  metric  for  actuator  force  output.  As 

mentioned  previously  the  performance  of  a  control  system  must 

be  implemented  with  real  actuators  that  have  a  limited 

ability  to  apply  control  energy.  A  metric  based  upon  either 

maximum  force  of  any  one  actuator  or,  possibly,  average  energy 

supplied  by  an  actuator  could  be  used  to  screen  modal 

arrangements  and  select  controller  gain.  This  second  metric 

should  be  used  in  conjunction  with  the  response  time  metric 

and  the  signal -to-noise  type  metric.  This  second  metric  may 

also  be  of  benefit  in  selecting  the  range  of  weighting 

matrices,  R  and  Q  used  in  calculating  gains.  The  important 
o  o 

decision  as  to  what  range  to  investigate  the  system  affects 
the  final  results  since  the  techniques  employed  search  for  a 
local  optimum  and  not  a  global  optimum. 

The  third  additional  recommendation  for  further  study  is 
to  investigate  other  sources  of  noise  in  the  design  process. 


The  effects  of  perturbing  the  structural  elements  should  be 
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examined.  A  change  to  the  stiffness  of  an  element  will 
change  the  location  of  the  natural  frequencies,  thus  inducing 
noise  to  the  system.  Calico  (3:44)  investigated  structural 
perturbations  on  the  CSDL-I  by  examining  open  loop  criteria 
in  the  form  of  reintroduction  of  coupling  within  the  control 
and  observation  matrices.  His  preferred  modal  arrangement, 
345-12-678,  was  robust  to  structural  perturbations,  but  not 
extremely  robust  to  sensor  noise.  This  would  indicate  that  a 
good  design  would  be  a  comprise  in  robustness  to  various 
types  of  noise  like  sensor  noise  and  structural 
miss-specification  noise.  This  investigation  could  easily  be 
accomplished  by  including  the  last  lines  of  code  from 
STIFF. M,  Appendix  B  into  the  beginning  of  DECOUPLE. M, 

Appendix  F  and  writing  additional  code  to  mesh  the  two 
together . 

The  final  additional  recommendation  for  further  study  is 
to  build  and  exercise  a  control  system  based  upon  the 
proposed  design  techniques.  The  true  test  of  a  good  design 
must  be  done  in  the  lab.  The  design  procedures  only  used 
knowledge  of  the  first  eight  modes.  The  four  residual  modes 
represent  all  uncontrolled,  unknown  modes.  It  is  asserted 
that  if  the  controller  doesn't  destabilize  the  four  residual 
modes  it  will  probably  not  destabilize  any  other  uncontrolled 
modes  present.  To  verify  this  and  other  principles  an  actual 
experiment  should  be  performed. 


Appendix  A 

Global  Coordinate  Martrices  K  and 
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Extended  Stiffness  Matrix  for  Element  3  in  s  Frame 

Columns 


1 

2 

3 

4 

5 

6 

7 

8 

9 

^11 

0.2500 

-0.1443 

-0.4082 

0 

0 

0 

-0.2500 

0.1443 

0.4082 

0 

-0.1443 

0.0833 

0.2357 

0 

0 

0 

0.1443 

-0.0833 

-0.2357 

-0.4082 

0.2357 

0.6667 

0 

0 

0 

0.4082 

-0.2357 

-0.6667 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.2500 

0.1443 

0.4082 

0 

0 

0 

0.2500 

-0.1443 

-0.4082 

0.1443 

-0.0833 

-0.2357 

0 

0 

0 

-0.1443 

0.0833 

0.2357 

0.4082 

-0.2357 

-0.6667 

0 

0 

0 

-0.4082 

0.2357 

0.6667 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Input  Location  Matrix,  D  in  s  Frame 

Col  ..ru  s 

1  2  3  5 _ 6 


0 

0 

0 

-0.3536 

0.6124 

-0.7071 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0.3536 

-0.6124 

-0.7071 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

-0.3536 

-0.6124 

-0.7071 

0 

O 

0 


0 

0 

0 

0 

0 

0 

0.3536 

0.6124 

-0.7071 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0.7071 

0 

-0.7071 


0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.7071 

0 

-0.7071 


Total  Stiffness  Matrix  K  in  s  Frame  with  Boundry  Conditions 

Columns 


1 

2 

3 

4 

5 

6 

13.7500 

6.4952 

18.3712 

-12.5000 

-7.2169 

-20.4124 

6.4952 

6.2500 

10.6066 

-7.2169 

-4.1667 

-11.7851 

18.3712 

10.6066 

40.0000 

-20.4124 

-11.7851 

-33.3333 

-12.5000 

-7.2169 

-20.4124 

79.4194 

21.2129 

20.4124 

-7.2169 

-4.1667 

-11.7851 

21.2129 

54.9249 

11.7851 

-20.4124 

-11.7851 

-33.3333 

20.4124 

11.7851 

51.0110 

-1.2500 

0.7217 

2.0412 

-50.0000 

0 

0 

0.7217 

-0.4167 

-1.1785 

0 

0 

0 

2.0412 

-1.1785 

-3.3333 

0 

0 

0 

0.0000 

0.0000 

0.0000 

-12.5000 

-21.6506 

0 

0.0000 

-1.6667 

2.3570 

-21.6506 

-37.5000 

0 

0.0000 

2.3570 

-3.3333 

0 

0 

0 

7 

8 

Columns 

9 

Continued 

10 

11 

12 

-1.2500 

0.7217 

2.0412 

0.0000 

0.0000 

0.0000 

0.7217 

-0.4167 

-1.1785 

0.0000 

-1.6667 

2.3570 

2.0412 

-1.1785 

-3.3333 

0.0000 

2.3570 

-3.3333 

-50.0000 

0 

0 

-12.5000 

-21.6506 

0 

0 

0 

0 

-21.6506 

-37.5000 

0 

0 

0 

0 

0 

0 

0 

68.1694 

-14.7177 

-2.0412 

-12.5000 

21.6506 

0 

-14.7177 

51.1749 

1.1785 

21.6506 

-37.5000 

0 

-2.0412 

1.1785 

21.0110 

0 

0 

0 

-12.5000 

21.6506 

0 

42.6777 

0.0000 

0.0000 

21.6506 

-37.5000 

0 

0.0000 

76.6667 

-2.3570 

0 

0 

0 

0.0000 

-2.3570 

21.0110 

Appendix  B 
Listing  of  STIFF. M 
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%THIS  IS  A  COMMAND  FILE  FOR  FINDING  THE  STIFFNESS  MATRIX 
% INITIAL  VARIABLES 

%ZR-ZERO  ROW , ZC-ZERO  COLUMN , L-ROTATION  MATRICES ,K- STIFFNESS 
Z  R  =  [  0 , 0,0, 0,0,0; 0,0,0, 0,0,0 ; 0,0,0, 0,0,0]; 

ZC=[ZR,ZR,ZR,ZR,ZR] ; 

ZC=ZC *  * 

K=[l,o[o,-l,0,0;ZR((l:2), : );-l,0,0,l,0,0;ZR( (1:2) , : )]; 

%FIRST  FIND  THE  INDIVIDUAL  ROTATION  MATRICES 
%DEFINE  MACROS 

R2= ’ [ cos (B) ,0,-sin(B) ;0,l,0;sin(B) ,0,cos(B) ) * ; 

R3= ’ [cos(C) ,sin(C) ,0; -sin(C) , cos(C) ,0;0 ,0 , 1] ' ; 

%FINDING  ROTATION  MATRIX  FOR  ELEMENT  ONE,  B  AND  C  ARE  ANGLES 
B=-. 95531661830; 

C=5*pi/6; 

Ll=eval (R3)*eval (R2) ; 

L1=[L1' ,ZR( : ,(1:3));ZR(: ,(1:3)), LI']; 

%FINDING  ROTATION  MATRIX  FOR  ELEMENT  TWO 
L2=[l, 0,0; 0,1,0; 0,0,1]; 

L2=[L2' ,ZR(: ,(1:3));ZR(: , (1:3) ) ,L2’ ] ; 

%FINDING  ROTATION  MATRIX  FOR  ELEMENT  THREE 
C=pi/6; 

L3=eval (R3)*eval (R2 ) ; 

L3=[L3',ZR(: ,(1:3));ZR(: , (1 : 3) ) , L3 ' ] ; 

%FINDING  REMAINING  ROTATION  MATRICES 
C=-pi/2; 

L4=eval (R3)*eval (R2) ; 

L4=[L4' ,ZR(: ,(l:3));ZR(: , (1:3) ) ,L4’ ) ; 

C=-pi/3 ; 

L5=eval (R3 ) ; 

L5=[L5' ,ZR(: ,(1:3));ZR(: , ( 1 : 3) ) ,L5 ' ] ; 

C=-2*pi/3; 

L6=eval (R3) ; 

L6= [L6 ' ,ZR( : , ( 1 : 3 ) ) ; ZR( : ,(1:3)),L6’]; 

B=-pi/4; 

L7=eval (R3)*eval (R2) ; 

L7=[L7' ,ZR(: ,(l:3));ZR(: , (1:3) ) ,L7* ] ; 

C=pi/3; 

L8=eval (R3)*eval (R2) ; 

L8=[L8' ,ZR(: ,(1:3));ZR(: , (1: 3) ) ,L8 ' ] ; 

C=2*pi/3 ; 

L9=eval (R3)*eval (R2 ) ; 

L9=[L9’ ,ZR(: ,(1:3));ZR(: , (1 : 3 ) ) ,L9 ’ ] ; 

C=-pi/3; 

L10=eval (R3)*eval (R2) ; 

L10= [L10  * ,ZR(: ,(1:3));ZR(: , (1: 3) ) ,L10 ’ ] ; 

Lll=eval (R2 ) ; 

Lll= [Lll * , ZR( : ,(1:3)) ;ZR( : , (1 : 3) ) , Lll ' ] ; 

C=pi ; 

Ll2=eval (R3)*eval (R2) ; 

L12=[L12' ,ZR(: ,(1:3));ZR(: , ( 1 : 3 ) ) , L12 ' ] ; 

%FIND  THE  ROTATED  STIFFNES  MATRICES  FOR  EACH  ELEMENT 
K1=L1'*K*L1; 

K2=L2 ' *K*L2 ; 
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K3=L3'*K*L3; 

K4=L4'*K*L4; 

K5=L5'*K*L5; 

K6=L6'*K*L6; 

K7=L7'*K*L7; 

K8=L8’*K*L8; 

K9=L9  *  *K*L9 ; 

K10=L10  *  *K*L10 ; 

Kll-Lll ' *K*L11 ; 

K12=L12'*K*L12; 

%FIND  EXTENDED  STIFFNESS  MATRICES  FOR  GLOBAL  COORDINATES 
K1E= [ Kl ; ZR ; ZR ; ZR ; ZR ; ZR ; ZR ; ZR ; ZR ] ; 

K1E= [KlE , ZC , ZC , ZC , ZC , ZC , ZC , ZC , ZC ] ; 

K2E= [ ZR ; K2 ; ZR ; ZR ; ZR ; ZR ; ZR ; ZR ; ZR ] ; 

K2E=[ZC,K2E,ZC,ZC,ZC,ZC , ZC , ZC , ZC] ; 

K3E=[K3( (1:3) , : ) ;ZR;K3( (4 : 6) , : ) ;ZR;ZR;ZR;ZR;ZR;ZR;ZR] ; 
K3E=[K3E( : , (1:3)), ZC , K3E( : , ( 4 : 6 ) ) , ZC, ZC, ZC , ZC, ZC, ZC , ZC] ; 
K4E= [K4( (1:3), : ) ;ZR;ZR;K4( (4: 6) , : ) ;ZR;ZR;ZR;ZR;ZR;ZR] ; 
K4E=[K4E( : , (1:3) ) ,ZC,ZC,K4E( : , (4: 6) ) ,ZC,ZC,ZC,ZC,ZC,ZC] ; 
K5E= [ ZR ; K5 ( ( 1 : 3 ) , : ) ; ZR;K5( (4 : 6) , : ) ;ZR;ZR; ZR; ZR; ZR ; ZR] ; 
K5E=[ZC,K5E( : , ( 1 : 3 ) ) , ZC ,K5E( : , ( 4 : 6) ) , ZC , ZC , ZC , ZC , ZC , ZC] ; 
K6E= [ ZR ; ZR ; K6 ; ZR ; ZR ; ZR ; ZR ; ZR ; ZR ] ; 

K6E=  [  ZC,  ZC  ,K6E,  ZC,  ZC  ,  ZC  ,  ZC,  ZC ,  ZC]  ; 

K7E= [ ZR ; K7 ( ( 1 : 3 ) , : ) ; ZR ; ZR ; K7 ( ( 4 : 6 ) , : ) ; ZR ; ZR ; ZR ; ZR ; ZR ] ; 
K7E=[ZC,K7E( : , (1 : 3 ) ) , ZC,ZC,K7E( : , ( 4 : 6) ) ,ZC ,ZC, ZC, ZC, ZC] ; 
K8E= [ ZR ; K8 ( ( 1 : 3 ) , : ) ;ZR;ZR;ZR;K8( (4: 6) , : ) ;ZR;ZR;ZR;ZR] ; 
K8E=[ZC,K8E( : , ( 1 : 3 ) ) , ZC, ZC, ZC ,K8E( : , (4: 6) ) ,ZC,ZC,ZC,ZC] ; 
K9E=[ZR;ZR;K9((l:3) , : )  ;ZR;ZR;ZR,*K9((4:6) , : )  ;ZR;ZR;ZR]  ; 

K9E= [ ZC , ZC , K9E( : , ( 1 : 3 ) ) , ZC , ZC , ZC , K9E( :,(4:6)),ZC,ZC,ZC]; 
K10E= [ ZR ; ZR ; K10 ((1:3), : ) ;ZR;ZR;ZR;7R;K10( (4:6) , :);ZR;ZR]; 
K10E= [ZC , ZC,K10E( : , ( 1 : 3) ) , ZC , ZC, ZC , ZC ,K10E( :,(4:6)),ZC,ZC]; 
K11E=[ZR;ZR;ZR;K11( (1:3) , : ) ;ZR;ZR;ZR;ZR;K11( (4: 6) , : ) ;ZR] ; 
K11E= [ ZC , ZC , ZC , K11E ( : , ( 1 : 3 ) ) , ZC , ZC , ZC , ZC ,K11E( : , ( 4 : 6) ) , ZC] ; 
K12E=[ZR;ZR;ZR;K12( (1:3), : ) ;ZR;ZR;ZR;ZR;ZR;K12( (4:6) ,:)]; 
K12E=[ZC,ZC,ZC,K12E( : , ( 1 : 3 ) ) , ZC , ZC , ZC , ZC, ZC ,K12E( : ,(4:6))]; 
%COMBINE  THE  STIFFNESS  MATRICES  TOGETHER , THEN  APPLY  BOUNDRY 
%COND.  THE  BOUNDRY  CONDITONS  ARE  U13-U30  =  0 
%SINCE  ONLY  INTERESTED  IN  DYNAMIC  RESPONSE  APPLY  BOUNDRY 
%COND  THEN  ADD  STIFFNESS  TOGETHER,  THIS  CREATES  THE  FORM 
%K  TOTAL =  EA/L*K1E  +  EA/L*K2E  ETC.  WHERE  EACH  EA/L  COULD  BE 
%DIFFERENT  FOR  EACH  ELEMENT,  APPLY  NOISE  TO  THE  EA/L  TERM 
K1D=K1E( (1:12) , (1:12)) ; 

K2D=K2E( ( 1 : 12 ) , ( 1 : 12 ) ) ; 

K3D=K3E( (1:12), (1:12)); 

K4D=K4E( (1:12), (1:12)); 

K5D=K5E( (1:12), (1:12)); 

K6D=K6E( (1:12) , (1:12)); 

K7D-K7E( (1:12), (1:12)); 

K8D=K8E( (1:12), (1:12)) ; 

K9D=K9E((1:12) ,(1:12)); 

K10D=K10E( ( 1 : 12 ) , ( 1 : 12 ) ) ; 

K11D=K11E( (1:12), (1:12)); 

K12D=K12E( (1:12), (1:12)); 

%COMBINE  TO  TOTAL  FINAL  STIFFNESS 
Cl=100 ; 
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C2=10; 

C3=100/sqrt(8) ; 

LEGS=C3* (K7D+K8D+K9D+K10D+K11D+K12D) ; 
KD=C1*K1D+C1*K2D+C2*K3D+C2*K4D+C1*K5D+C1*K6D+LEGS 


Appendix  C 
Listing  of  SENSOR. M 
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%THIS  IS  A  COMMAND  FILE  FOR  FINDING  THE  SENSOR  MATRIX  D 
%INITIAL  VARIABLES 

%ZR-ZERO  ROW , ZC-ZERO  COLUMN , L -ROTATION  MATRICES 
ZR= [0,0,0,0,0,0;0,0,0,0,0,0;0, 0,0, 0,0,0]; 
ZC=[ZR,ZR,ZR,ZR,ZR] ; 

ZC=ZC'  ; 

K=[1,0,0,  0 ,0 ,0;ZR( (1: 2) , : ) ;  0,0 ,0 ,0 ,0 ,0;ZR( (1:2) , : ) ] ; 
%FIRST  FIND  THE  INDIVIDUAL  ROTATION  MATRICES 
%DEFINE  MACROS 

R2= ' [ cos(B) ,0,-sin(B);0,l,0;sin(B),0, cos (B) ] ' ; 

R3= ' [cos(C) ,sin(C),0;-sin(C) ,cos(C) ,0;0,0,1] ' ; 

C=-2*pi/3; 

B=-pi/4; 

L7=eval (R3)*eval (R2 ) ; 

L7=[L7' ,ZR(: ,(1:3));ZR(: ,(1:3)),L7' ]; 

C=pi/3; 

L8=eval (R3)*eval (R2) ; 

L8=tL8’ ,ZR(: ,(1:3));ZR(: , (1 : 3) ) ,L8  '  ] ; 

C=2*pi/3; 

L9=eval (R3)*eval (R2) ; 

L9=[L9‘ ,ZR( : ,(l:3));ZR(: , (1 : 3) ) ,L9’ ] ; 

C=-pi/3; 

L10=eval(R3)*eval(R2) ; 

L10=[L10' ,ZR( : , (1:3)); ZR( : , (1 : 3) ) ,L10 ’ ] ; 

Lll=eval (R2 ) ; 

Lll= [Lll ' ,ZR(: ,(l:3));ZR(: , (1 : 3) ) ,L11 ’ ) ; 

C=pi; 

L12=eval(R3)*eval(R2); 

L12=[L12’ ,ZR(: , (l:3));ZR(: , ( 1 : 3) ) , L12 * ] ; 

K7=L7’*K; 

K8=L8 ' *K; 

K9=L9'*K; 

K10=L10**K; 

K11=L11’*K; 

K12=L12 ' *K; 

K7E= [ZR;K7 ( (1 : 3 ) , : ) ; ZR; ZR;K7 ( ( 4 : 6) , : ) ; ZR; ZR; ZR;ZR;ZR] ; 
K7E=[ZC,K7E( : , (1 : 3) ) ,ZC, ZC,K7E( : , ( 4 : 6) ) ,ZC,ZC,ZC,ZC,ZC] ; 
K8E=[ZR;K8( (1:3) , : ) ; ZR; ZR; ZR; K8 ( ( 4 : 6) , : ) ;ZR;ZR;ZR;ZR] ; 

K8E= (ZC,K8E( : , ( 1 : 3 ) ) , ZC , ZC , ZC , K8E( : , ( 4 : 6) ) , ZC , ZC , ZC , ZC] ; 
K9E=[ZR;ZR;K9((1:3), : ) ;ZR;ZR;ZR;K9( (4: 6) , : ) ;ZR;ZR;ZR] ; 

K9E= [ZC , ZC , K9E( : , (1 : 3 ) ) , ZC , ZC, ZC ,K9E( :,(4:6)),ZC,ZC,ZC); 
K10E= [ ZR ; ZR ; K10 ( ( 1 : 3 ) , : ) ; ZR ; ZR ; ZR ; Z? *  K10 ((4:6),:) ; ZR ; ZR ] ; 
K10E= [ZC , ZC ,K10E( : , ( 1 : 3) ) , ZC , ZC , ZC , ZC ,K10E( : , ( 4: 6) ) ,ZC,ZC] ; 
K11E= [ ZR ; ZR ; ZR ; Kll ( ( 1 : 3 ) , : ) ; ZR ; ZR ; ZR ; ZR ; Kll ( ( 4 : 6 ) , : ) ; ZR ] ; 
K11E= [ ZC , ZC , ZC , K11E( : , ( 1 : 3 ) ) , ZC , ZC , ZC , ZC ,K11E( : , ( 4 : 6) ) , ZC] ; 
K12E=[ZR;ZR;ZR;K12( (1 : 3) , : ) ;ZR;ZR;ZR;ZR;ZR;K12( (4: 6) , : )]; 
K12E= [ZC , ZC , ZC ,K12E( : , (1: 3) ) ,ZC,ZC,ZC,ZC,ZC,K12E( : ,(4:6))]; 
K7D=K7E( (1 : 12) , ( 1 : 12 ) ) ; 

K8D=K8E( (1:12), (1:12)); 

K9D=K9E( (1:12), (1:12)); 

K10D=K10E( (1:12), (1:12)); 

K11D=K11E( (1:12) ,(1:12)); 

K12D=K12E( (1:12), (1:12)); 
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%COMB INI MG  INTO  PINAL  MATRIX 

D= [K"D( : ,4) ,K8D( : ,4) ,K9D( : ,7) ,K10D( : ,7) , 
K11D( : ,10) ,K12D( : ,10)]; 


State-space  Matrices  A,B  and  C 


State-space  Plant  Matrix  A  Full-order  System 

Columns 

2  3  4  5  6 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

•1.8010 


0  -2.7714 


0  -8.3562 


0  -8.7462 


0  -11.5474 


0  -17.6777 


0  -21.7349 


D-2 


Columns  Continued 


8 

9 

10 

11 

12 

13 

14 

0 

0 

0 

0 

0 

1. 

0000 

0 

0 

0 

0 

0 

0 

0 

1. 

0000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0. 

0134 

0 

0 

0 

0 

0 

0 

0 

-0.1 

0166 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-22.6127 

0 

0 

0 

0 

0 

0 

0  -72. 

9227 

0 

0 

0 

0 

0 

0 

0  - 

85.5732 

0 

0 

0 

0 

0 

0 

0 

-105.7785 

0 

0 

0 

0 

0 

0 

0  -166.5441 

0 

0 

Columns  Continued 

15 

16 

17 

18  19 

20 

21 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1.0000 

0 

0 

0 

0 

0 

0 

0  1 

.0000 

0 

0 

0 

0 

0 

0 

0 

1.0000 

0 

0 

0 

0 

0 

0 

0 

1.0000 

0 

0 

0 

0 

0 

0 

0  1.0000 

0 

0 

0 

0 

0 

0 

0  1.0000 

0 

0 

0 

0 

0 

0 

0 

1.0000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.0289 

0 

0 

0 

0 

0 

0 

0  -0 

.0296 

0 

0 

0 

0 

0 

0 

0 

-0.0340 

0 

0 

0 

0 

0 

0 

0 

-0.0420 

0 

0 

0 

0 

0 

0 

0  -0.0466 

0 

0 

0 

0 

0 

0 

0  -0.0476 

0 

0 

0 

0 

0 

0 

0 

-0.0854 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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Columns  Continued 


22  23  24 _ 

0  0  0 

0  0  0 

0  0  0 

0  0  0 

0  0  0 

0  0  0 

0  0  0 

0  0  0 

0  0  0 

1.0000  0  0 

0  1.0000  0 

0  0  1.0000 

0  0  0 

0  0  0 

0  0  0 

0  0  0 

0  0  0 

0  0  0 

0  0  0 

0  0  0 

0  0  0 

-0.0925  0  0 

0-0.1028  0 
0  0-0.1291 


State-space  Input  Matrix  B 
Columns 

1  2 _ 3 _ 4 _ 5  6 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0621 

-0.0976 

-0.0651 

-0.3518 

0.4969 

-0.4082 

0.0689 

0.0978 

0.3266 

0.4483 

-0.3117 

-0.1611 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.0621 

-0.0976 

-0.0651 

0.3518 

0.4969 

0.4082 

-0.0689 

0.0978 

0.3266 

-0.4483 

-0.3117 

-0.1611 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.0948 

-0.0236 

-0.3838 

0.0843 

-0.0686 

-0.4082 

-0.5217 

-0.4232 

0.3539 

-0.2123 

0.2067 

0.0187 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.0327 

0.1582 

0.1095 

-0.2674 

0.2200 

0.4082 

-0.4528 

-0.5162 

-0.3238 

0.2360 

-0.2049 

-0.0351 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0327 

0.1582 

0.1095 

0.2674 

0.2200 

-0.4082 

0.4528 

-0.5162 

-0.3238 

-0.2360 

-0.2049 

-0.0351 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0948 

-0.0236 

-0.3838 

-0.0843 

-0.0686 

0.4082 

0.5217 

-0.4232 

0.3539 

0.2123 

0.2067 

0.0187 


State-space  Output  Matrix,  C 
Columns 


0.0621  -0.0976  -0.0651  -0.3518  0.4969  -0.4082  0.0689 

-0.0621  -0.0976  -0.0651  0.3518  0.4969  0.4082  -0.0689 

-0.0948  -0.0236  -0.3838  0.0843  -0.0686  -0.4082  -0.5217 

-0.0327  0.1582  0.1095  -0.2674  0.2200  0.4082  -0.4528 

0.0327  0.1582  0.1095  0.2674  0.2200  -0.4082  0.4528 

0.0948  -0.0236  -0.3838  -0.0843  -0.0686  0.4082  0.5217 


Columns  Continued 

S  9  10  11  12  13-24 


0.0978 

0.0978 

-0.4232 

-0.5162 

-0.5162 

-0.4232 


0.3266 

0.3266 

0.3539 

-0.3238 

-0.3238 

0.3539 


0.4483 

-0.4483 

-0.2123 

0.2360 

-0.2360 

0.2123 


-0.3117 

-0.3117 

0.2067 

-0.2049 

-0.2049 

0.2067 


-0.1611 

-0.1611 

0.0187 

-0.0351 

-0.0351 

0.0187 


0.  .  .0 


0. 


.0 
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I 

► 

► 
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Total  Plant  Matrix  of  Decoupled  Control 


State-space  Closed  Looped  Plant  Matrix,  Ae,  Suppressed 

for  Arrangement  123-45-678  Q  =20,  R  =1 

o  o 

Columns 


1 

2 

3 

4 

5 

6 

7 

0 

0 

0 

1.0000 

0 

0 

0 

0 

0 

0 

0 

1.0000 

0 

0 

0 

0 

0 

0 

0 

1.0000 

0 

-1.9495 

0.0000 

0.0000 

-0.9242 

0.0000 

0.0000 

-0.1485 

0.0000 

-2.8973 

-0.4290 

0.0000 

-1.3451 

-0.4175 

0.0000 

0.0000 

0.3469 

-9.0286 

0.0000 

-0.7197 

-2.6710 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0 

0 

0 

0 

0 

0 

-1.8010 

0 

0 

0 

0 

0 

0 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.2250 

0.0000 

0.0000 

1.3797 

0.0000 

0.0000 

0.2250 

0.0000 

0.1455 

0.0873 

0.0000 

0.5705 

-0.3528 

0.0000 

-0.4261 

0.0000 

0.0000 

0 

0 

0 

0 

0.0000 

-0.2210 

0.3103 

0 

0 

0 

0 

-0.0513 

0.0000 

0.0000 

0 

0 

0 

0 

0.0000 

-0.0238 

0.0334 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.7314 

0.0000 

0.0000 

-4.4856 

0.0000 

0.0000 

-0.7314 

0.0000 

0.9032 

0.6338 

0.0000 

3.7579 

-2.0031 

0.0000 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0.6323 

0.0000 

0.0000 

0 

0 

0 

0 

0.0000 

-0.7694 

0.9432 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0.0331 

0.0000 

0.0000 

0 

0 

0 

0 

0.0000 

-0.0348 

0.0426 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0000 

-0.0464 

1.3373 

0.0000 

3.0425 

2.9039 

0.0000 

-0.4294 

0.0000 

0.0000 

-2.6332 

0.0000 

0.0000 

-0.4294 

0.0000 

-0.2373 

0.2378 

0.0000 

-0.0323 

1.3528 

0.0000 

0.0000 

-0.0684 

-0.0992 

0.0000 

-0.4055 

0.0468 

0.0000 

E-2 


8 

9 

10 

Columns 

11  12 

13 

14 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

-0.9108 

0.0000 

0.0000 

0.0000 

0.0000 

-0.1259 

-0.4290 

0.0000 

-1.3284 

-0.4175 

0.0000 

0.0000 

0.3469 

-0.6724 

0.0000 

-0.7197 

-2.6421 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0 

0 

0.0000 

0.0000 

-0.9230 

-2.4710 

0 

1.0000 

0 

0.0000 

0.0000 

-0.8526 

-2.2825 

0 

0 

1.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0134 

0 

0 

0.0000 

0.0000 

-2.1047 

1.7849 

0 

-0.0166 

0 

0.0000 

0.0000 

-0.9802- 

10.9804 

0 

0 

-0.0289 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

1.37  97 

0.0000 

0.0000 

-9.1275 

0.0000 

0.1455 

0.0873 

0.0000 

0.5705 

-0.3528 

0.0000- 

-12.0382 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-4.4856 

0.0000 

0.0000 

-0.3508 

0.0000 

0.9032 

0.6338 

0.0000 

3.7579 

-2.0031 

0.0000 

0.1737 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0.4875 

0.0000 

0 

0 

0 

0 

0 

0.0000 

-0.3407 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0.0255 

0.0000 

0 

0 

0 

0 

0 

0.0000 

-0.0154 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.0464 

1.3373 

0.0000 

3.0425 

2.9039 

0.0000 

-0.1001 

0.0000 

0.0000 

-2.6332 

0.0000 

0.0000 

0.3978 

0.0000 

-0.2373 

0.2378 

0.0000 

-0.0323 

1.3528 

0.0000 

0.3391 

-0.0684 

-0.0992 

0.0000 

-0.4055 

0.0468 

0.0000 

0.1418 

E-3 


Columns 


15 

16 

17 

18 

19 

20 

21 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

1.0000 

0 

0 

0 

0 

0 

0 

0 

1.0000 

0 

0 

0 

0 

0 

-2.7531 

0.0000 

-0.3814 

0.0000 

-2.7236 

0.0000 

0.0000 

0.0000 

-3.5440 

0.0000 

-0.4909 

0.0000 

-3.5101 

0.0000 

0 

0 

-2.9269 

0.0000 

1.0000 

0 

0.0000 

0 

0 

0.0000 

-3.5573 

0 

1.0000 

0.0000 

0 

0 

-9.0985 

0.0000 

-0.0296 

0 

0.0000 

0 

0 

0.0000- 

11.9308 

0 

-0.0340 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000- 

18.2346 

-2.5054 

0.0000 

-0.3508 

0.0000 

-2.5054 

0.0000 

0.0000 

0.0000 

1.2424 

0.0000 

0.1737 

0.0000 

1.2424 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0000 

-0.7155 

0.0000 

-0.1001 

0.0000 

-0.7155 

0.0000 

2.8408 

0.0000 

0.3978 

0.0000 

2.8408 

0.0000 

0.0000 

0.0000 

2.4249 

0.0000 

0.3391 

0.0000 

2.4249 

0.0000 

0.0000 

1.0139 

0.0000 

0.1418 

0.0000 

1.0139 

0.0000 

E-4 


Columns 


22 

23 

24 

25 

26 

27 

28 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

1.0000 

0 

0 

0 

0 

0 

0 

0 

1.0000 

0 

0 

0 

0 

0 

0 

0 

1.0000 

0 

0 

0.0000 

0.0000 

-4.5952 

0.0000 

0.0000 

-0.5569 

0.0000 

21.7349 

0.0000 

0.0000 

-0.0466 

0.0000 

0.0000 

0.0000 

0.0000- 

22.6127 

0.0000 

0.0000 

-0.0476 

0.0000 

0.0000 

0 

0 

0 

0 

0 

-4.5533 

0.0000 

0 

0 

0 

0 

0 

0.0000 

-4.4433 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0 

0 

0 

0 

0- 

18.0441 

0.0000 

0 

0 

0 

0 

0 

0.0000- 

21.9678 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Columns 


29 

30 

31 

32 

33 

34 

35 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

-3.7106 

0.0000 

-1.0634 

0 

0 

0 

0 

-3.4276 

0.0000 

-0.9823 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

2.6804 

0.0000 

0.7682 

0 

0 

0 

0 

-3.9408 

0.0000 

-1.1294 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0.0000 

-3.6116 

0.0000 

0 

0 

0 

0 

0.8460 

0.0000 

-2.5745 

0 

0 

0 

0 

0.0000 

-0.4347 

0.0000 

0 

0 

0 

0 

0.0912 

0.0000 

-0.2775 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0000 

-4.5531 

0.0000 

0.0000 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

1.0000 

0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0 

1.0000 

0 

0.0000 

0.3205 

0.0000 

-4.3123 

0 

0 

1.0000 

0.4677 

0.0000 

-0.1155 

0.0000 

-0.0420 

0 

0 
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Listing  of  DECOUPLE.** 


P-1 


%THIS  M-PILE  CALCULATES  A  TIME  HISTORY  RESPONSE  METRIC  FOR  3 
%DECOUPLED  LQR  CONTROLLERS  APPLIED  TO  THE  CSDL1 
%THE  VARIABLES  REQUIRED  IN  THE  WORK  SPACE  ARE: 

%THE  IMPORTANT  OUTPUT  VARIABLES  ARE: 
toad  input. mat 

chdir  ga90m: [ jbrewer .matlab. output] 
for  i=l : length(n(l , : ) ) 

CONlMl=n(l,i) ; 

CONlM2=n(2 ,  i ) ; 

CONlM3=n( 3 ,  i ) ; 

CON2Ml=n( 4 ,  i  )  ; 

CON2M2=n(5,i) ; 

CON3Ml=n( 6 ,  i  ) ; 

CON3M2=n(7 ,i) ; 

CON3M3=n(8,i); 

RESMl=n( 9, i ) ; 

RESM2=n( 10 , i ) ; 

RESM3=n ( 1 1 ,  i ) ; 

RESM4=n(12,i) ; 

%FORMING  THE  INDIVIDUAL  A  MATRICES  FOR  EACH  CONTROLLER, 

%A  MODAL  DAMPING  OF  .005  IS  USED 

W=[V(CONlMl,l) ,0,0;0,V(CON1M2, 1) , 0 ; 0 , 0 , V(C0N1M3 , 1 ) ] ; 

ZW=  -2*  .  005*sqrt  (W)  ; 

W=-W  ; 

Al=[zeros(3) ,eye(3) ;W,ZW] ; 

W=  [V(CON2Ml ,  1 )  ,0  ,*0  ,  V(CON2M2 ,  1)  ]  ,* 

ZW=-2*.005*sqrt(W) ; 

W=-W; 

A2=[zeros(2) , eye(2) ;W,ZW] ; 

W=[V(CCN3M1,1),0,0;0,V(CON3M2,1),0;0,0,V(CON3M3,1)]; 
ZW=-2*.005*aqrt(W) ; 

W=-W ; 

A3=[aeros(3) ,eye(3) ;W,ZW]; 

%FIND  THE  RESIDUAL  A  MATRIX 

W=[V(RESM1,1) ,0,0,0;0,V(RESM2,1) , 0 , 0 ; 0 , 0 , V(RESM3 , 1 ) ,0; . . 

0 , 0 , 0 , V ( RESM4 , 1 ) ] ; 

ZW=-2*.005*sqrt(W); 

W=-W; 

AR=[zeros(4) ,eye(4) ;W,ZW] ; 

%FIND  THE  INDIVIDUAL  CONTROLLER'S  Bj  MATRIX  USING  THE  ORG 
%STATESPACE  MODEL  D  MATRIX  (INPUT  MATRIX' 

B=U'*D; 

Bl=[zeros(3,6) ;B(C0N1M1, : ) ;B(CONlM2, : ) ;B(CONlM3 , : ) ] ; 

B2= [ zeros (2,6) ; B( CON2M1 , : ) ;B(CON2M2 , : ) ] ; 

B3=[zeros(3, 6) ;B(CON3Ml , : ) ;B(CON3M2 , : ) ;B(CON3M3 , : ) ] ; 

%FIND  THE  RESDUAL  MODE'S  B  MATRIX 

BR=[zeros(4,6);B(RESMl, : );B(RESM2, : );B(RESM3, : ) ;B(RESM4 , : ) ] ; 
%DIVIDE  THE  UNPRETURBED  C  MATRIX  INTO  CORRECT  Cj  MATRIX 
%FOR  EACH  CONTROLLER 
C=(U'*D)’ ; 

C1-[C( : , CON1M1 ) ,C( : ,CONlM2) ,C( : ,CONlM3)]; 

Cl* [Cl, zeros (6, 3)] ; 

C2-[C( : , CON2M1 ) ,C( : ,CON2M2)]; 


F-2 


C2=[C2 , zeros ( 6 ,  2 )  ]  ; 

C3=[C( : , CON3M1 ) ,C( : ,CON3M2) ,C( : ,CON3M3)] ; 

C3=[C3f zeros (6, 3) j ; 

CR=[C( : , RESM1 )  ,C( : ,RESM2) ,C( : ,RESM3) ,C( : ,RESM4)]; 

CR= [ CR , zeros (6,4)]; 

%CONVERT  INDIVIDUAL  CONTROLLER  MATRICES  TO  A  DECOUPLED  SUB- 
%SPACE, USING  THE  RELATIONSHIP  BiTj=0  AND  GAMiCj=0  i=l,2  j=2,3 
Tl=eye( 6, 6) ; 

T2=null (Bl) ; 

T3=null ( [B1;B2] ) ; 

GAMl=nul 1 ( [C2 , C3 ] ' ) ; 

GAM2=null(C3' ); 

GAM3=eye(6,6) ; 

B1STAR=B1*T1 ; 

B2STAR=B2*T2; 

B3STAR=B3*T3; 

C1STAR=GAM1**C1; 

C2STAR=GAM2'*C2; 

C3STAR=GAM3’*C3; 

%DESIGN  A  LINEAR  QUADRATIC  REGULATOR  AND  ESTIMATOR  FOR  EACH 
%CONTROLLER,  Rc-CONROL  WEIGHTING,  QcL  AND  QcL  STATE  WT. 

%NOTE  THE  SIX  DIFFERENT  SIZES  OF  THE  CONTROL  WEIGHTING 
%GAINS.  Rcl  IS  USED  TO  FIND  CONI  CONTROLLER  GAINS,  WHILE 
%Rf 1  IS  USED  TO  FIND  CONI  ESTIMATOR  GAINS 
eval ( [ 'COMB* , int2str(i ) , ' = zeros ( length(Qc(l , 

1 ength(NOISE( 1 ,:)));  ’]) 
for  qc=l : length(Qc( 1 , : ) ) 

Rcl*eye( 1 ength(Tl ’ *T1 ) ) ; 

Rc2=eye(Iength(T2’*T2)); 

Rc3=eye( length(T3 ' *T3) ) ; 

Rf l=eye( length(GAMl ' *GAM1) ) ; 

Rf 2=eye( length(GAM2 ' *GAM2) ) ; 

Rf 3=eye( length(GAM3 ' *GAM3) ) ; 

%ASSUME  NOISE  IS  DECOUPLED  IN  THE  MODAL  STATES 
GAML=eye( 6 , 6) ; 

GAMS=eye( 4 , 4) ; 

%STATE  CONTROLL  WEIGHTING  VARIATION 

%eval ( [ 'COMB* ,int2str(i) , 1 =  zeros ( length(Qc(l , 

length(NOISE(l ,:)));  ’]) 

%for  qc=l : length(Qc(l , : ) ) 

Qcfl=20*eye(6,6) ; 

Qcf2=20*eye(4,4); 

Qcf 3=20*eye( 6 , 6) ; 

%Qcl=diag(Qc( (1 : 3) ,qc) ,0) ; 

Qcl=[20*eye(3, 3) , zeros (3 , 3) ; zeros ( 3 ,3) , . . 
diag( [20 ,Qc(l ,qc) ,35] ) j ; 

%QcI=Qc(qc)*eye(6,6) ; 

%Qc 2=diag(Qc( (4:5),qc),0); 

Qc2=[35*eye(2,2) , zeros (2, 2) ; zeros (2, 2) , 35*eye(2 ,2) ] ; 
%Qc2=Qc(qc)*eye(4,4) ; 

%Qc3=diag(Qc( (4: 6) ,qc) ,0) ; 

Qc3=[35*eye(3,3) , zeros ( 3 , 3 ) ; zeros ( 3 , 3) , . . 
diag( [Qc( 2 ,qc) ,35,Qc(3,qc)])]; 

%Qc3=Qc(qc)*eye(6,6) ; 

%FIND  GAINS  FOR  CONTROLLER  AND  ESTIMATOR  USING  Iqr,  lge 
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%A  MEG  SIGN  IS  REQUIRED  ON  CONTROLLER  GAIN  DUE  TO  THE  WAY 
%MATLAB  DEFINES  FEEDBACK,  AND  THE  WAY  THE  THESIS  EQUATIONS 
%DEFINE  IT. ALSO  NOTE  GAINS  ARE  CALCULATED  FOR  UNPRETURBED 
% STATES 

Kcl=-lqr ( A1 , BISTAR ,Qcl ,Rcl ) ; 

Kf l=lqe(Al ,GAML,ClSTAR,Qcf 1 ,Rf 1) ; 

Kc2=-lqr(A2 , B2STAR,Qc2 , Rc2) ; 

Kf 2= lqe( A2 ,GAMS ,C2STAR ,Qc£2 , Rf 2 ) ; 

Kc3=-lqr(A3,B3STAR,Qc3,Rc3); 

Kf 3=lqe(A3 ,GAML ,C3STAR,Qcf 3 ,Rf 3) ; 

%USE  ORGINAL  INPUT  MATRIX  B  WITH  ADDTI IONAL  ERROR  STATES  ADD 
%AS  THE  REFERENCE  INPUT  MATRIX  TO  CONVERT  AN  IMPULSE  AT  ONE 
%OF  THE  LEGS  INTO  THE  MODAL  COORDINATE  SYSTEM 
Btotal=[B; zeros (28, 6) ] ; 

% 

%THE  ABOVE  Btotal  WAS  THE  REFERENCE  INPUT  MATRIX  USED  TO 
%GENERATE  ALL  OF  THE  THESIS  DATA.  IT  ACTUALLY  APPLIES  INITIAL 
%VELOCITY  CONDITIONS  INSTEAD  OF  DISPLACEMENT  CONDITIONS.  THE 
%FINAL  SELECTED  COMBINATION  CHANGES  VERY  LITTLE  SINCE  ALL 
%COMBI NATIONS  ARE  COMPARED  TO  THE  SAME  REF.  INPUT.  IF  FURTHER 
%WORK  IS  DONE  I  WOULD  SUGGEST  CHANGING  THE  ABOVE  LINE  TO 
%  [zeros (12 ,6) ;B; zeros(16, 6) ] ;  AS  THIS  MATCHES  EQ  3.19 
%ALSO  CHANGING  THE  MOVING  AVERAGE  CUT  OF  FROM  .05  TO  .01  IS 
%NECESSARY  SINCE  THE  NEW  Btotal  INPUTS  LESS  ENERGY 
% 

%EXTRACT  NOISE  VARIABLES  FROM  INPUT  MATRIX  NOISE 
% THREE  LEVELS  OF  NOISE  ARE  .8  1.0  AND  1.2  MULTIPLIERS 
%NOISE7  IS  NOISE  APPLIED  TO  THE  SENSOR  ON  LEG  7,  ETC 
for  j=l: length(NOISE(l , : ) ) 

NOISE7=NOISE(l, j); 

NOISE8=NOISE(2,j); 

NOISE9=NOISE( 3 , j) ; 

NOISE10=NOISE( 4 , j) ; 

NOISEll=NOISE(5, j); 

NOISE12=NOISE( 6, j) ; 

%THIS  FINDS  THE  C  MATRIX  AFTER  APPLYING  THE  NOISE  TO  THE 
%SENSORS  IN  LOCAL  COORDINATE  FRAME  AND  ROTATING  INTO  MODAL  CD 
%ZR-ZERO  ROW , ZC-ZERO  COLUMN , L -ROTATION  MATRICES ,K- INPUT  LOCAL 
ZR= [0,0,0,0,0,0;0, 0,0, 0,0,0; 0,0, 0,0, 0,0]; 

ZC=  [ZR  ,  ZR ,  ZR ,  ZR  ,  ZR]  ,* 

ZC=ZC’ ; 

K=[1,0,0,  0,0,0;ZR((1:2),:);  0 ,0 ,0 , 0 ,0 ,0 ;ZR( (1 : 2) , : )  ] ; 

%FIRST  FIND  THE  INDIVIDUAL  ROTATION  MATRICES 
%DEFINE  MACROS 

R2= ' [cos(Ba) ,0, -sin(Ba) ;0 , 1 , 0 ; sin(Ba) ,0 ,cos(Ba) ] ' ; 

R3= ' [cos(Ca) , sin(Ca) ,0;-sin(Ca) , cos(Ca) ,0;0,0,1] ' ; 

Ca=-2*pi/3; 

Ba=-pi/4; 

L7=eval (R3)*eval (R2) ; 

L7=[L7 ' ,ZR( : ,(1:3));ZR( : ,(1:3)),L7' ]; 

Ca=pi/3; 

L8=eval (R3)*eval (R2) ; 

L8= [L8 ’ ,ZR( : , ( 1 : 3 ) ) ; ZR( ; , ( 1 : 3 ) ) , L8 ' ]; 

Ca=2*pi/3; 

L9=eval (R3)*eval (R2) ; 
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L9= [L9 ' ,ZR(: , ( 1 : 3) ) ; ZR( : ,(1:3)), L9’); 

Ca=-pi/3; 

L10=eval (R3)*eval (R2) ; 

L10= [L10  * , ZR(:,(1:3));ZR(:,(1:3)) ,L10 ' ] ; 

Lll=eval (R2 ) ; 

L11S[L11 ' ,  ZR( : , (1 : 3) ) ;ZR( : ,(1:3)) ,L11' ] ; 

Ca=pi ; 

L12=eval (R3)*eval (R2) ; 

L12= [L12 * ,ZR(: ,(1:3));ZR(: , (1:3) ) , L12 ' ] ; 

K7=L7’*(NOISE7*K); 

K8=L8'*(NOISE8*K); 

K9=L9 1  * (NOISE9*K) ; 

K10=L10 ' *(NOISE10*K) ; 

K11=L11'*(N0ISE11*K); 

K12=L12 ' * (NOISE12*K) ; 

K7E= [ ZR ; K7 ( ( 1 : 3 ) , : ) ;ZR;ZR;K7( (4: 6) , : ) ;ZR;ZR;ZR;ZR;ZR] ; 
K7E=[ZC,K7E( : , (1 : 3) ) ,ZC,ZC,K7E( : , ( 4 : 6) ) , ZC, ZC,ZC, ZC, ZC] ; 

K8E= [ ZR ; K8 ( ( 1 : 3 ) , : ) ;ZR;ZR;ZR;K8( (4: 6) , : ) ;ZR;ZR;ZR;ZR] ; 

K8E= [ZC , K8E( : , ( 1 : 3 ) ) , ZC , ZC , ZC , K8E( : , ( 4 : 6 ) ) , ZC , ZC , ZC , ZC] ; 

K9E= [ ZR ; ZR ; K9 ( (1 : 3 ) , : ) ; ZR; ZR; ZR;K9( ( 4 : 6) , : );ZR;ZR;ZR] ; 

K9E= [ ZC , ZC,K9E( : , ( 1 : 3 ) > , ZC , ZC , ZC ,K9E( : , (4: 6) ) ,ZC,ZC,ZC] ; 
K10E=[ZR;ZR,*K10((1:3) , : ) ; ZR ; ZR ; ZR ; ZR ; K10 ((4:6), : )  ;ZR;ZR]; 
K10E=[ZC,ZC,K10E( : ,(1:3)) ,ZC, ZC, ZC, ZC,K10E( :,(4:6)),ZC,ZC]; 

K11E= [ZR ; ZR ; ZR;K11( (1 : 3) , : ) ;ZR;ZR;ZR;ZR;K11( (4: 6) , : ) ;ZR] ; 
K11E=[ZC,ZC,ZC,K11E( : , (1:3) ) ,ZC,ZC,ZC,ZC,K11E( : , (4: 6) ) f ZC] ; 
K12E=[ZR,*ZR;ZR;K12((l:3) , : )  ;ZR,*ZR;ZR;ZR;ZR;K12(  (4: 6)  ,:)]; 

K12E= [ ZC , ZC , ZC ,K12E( : ,(1:3)), ZC, ZC , ZC, ZC , ZC ,K12E( : ,(4:6))]; 
K7D=K7E( (1: 12) ,(1:12)); 

K8D=K8E( (1:12), (1:12)); 

K9D=K9E( (1:12) ,(1:12) ) ; 

K10D=K10E( (1:12), (1:12)); 

K11D=K11E( (1:12), (1:12)); 

K12D=K12E( (1:12) , (1:12) ) ; 

Dn=[K7D(: ,4),K8D( : ,4),K9D( : ,7),K10D( : ,7),K11D( : ,10), . . 

K12D( : , 10 ) ] ; 

Ctotal=[(U'*Dn) * , zeros (6, 28) ] ; 

%DIVIDE  THE  PRETURBED  C  MATRIX  INTO  CORRECT  Cj  MATRIX 
%FOR  EACH  CONTROLLER 

Cel=[Ctotal ( : , C0N1M1 ) ,Ctotal ( : ,C0N1M2 ) ,Ctotal ( : ,C0N1M3)]; 

Cel =[ Cel , zeros (6, 3) ] ; 

Ce2= [Ctotal ( : ,C0N2M1 ) ,Ctotal ( : ,CON2M2)]; 

Ce2=[Ce2 , zeros (6, 2) ] ; 

Ce3= [Ctotal ( : ,C0N3M1) ,Ctotal ( : ,CON3M2 ) ,Ctotal ( : ,CON3M3)]; 
Ce3*[Ce3,zeros(6,3) ] ; 

CeR= [Ctotal ( : , RESM1 ) , Ctotal ( : , RESM2 ), Ctotal (: , RESM3 ) , . . 

Ctotal (: ,RESM4) ] ; 

CeR= [CeR , zeros (6,4)]; 

4COMBINING  THE  DECOUPLED  CONTROLLERS  BACK  INTO  A  SINGLEPLANT 
%MATRIX  INCLUDING  THE  ERROR  STATES  (MODES  NOT  IN  ORDER  !) 

%NOTE  PRETURBED  C  MATRICES  ARE  USED  IN  FEEDBACK! 
Atotal=[(Al+BlSTAR*Kcl) , BlSTAR*Kcl ,Bl*T2*Kc2 ,Bl*T2*Kc2 , . . 

Bl*T3*Kc3 ,Bl*T3*Kc3 , zeros( 6 , 8) ; . . 
zeros (6, 6) , (Al-Kf 1*GAM1 ' *Cel ) ,Kfl*GAMl’*Ce2 ,zeros(6, 4) , . . 

Kf 1*GAM1 * *Ce3 , zeros (6,6) ,Kf 1*GAM1 ' *CeR] ; 

Atotal= [Atotal ;B2*Tl*Kcl ,B2*Tl*Kcl , (A2+B2STAR*Kc2) ,B2STAR*Kc2 , . . 
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B2*T3*Kc3 , B2*T3*Kc3 , zeros (4,8); . . 

Kf 2*GAM2 ’ *Cel , zeros (4,6), zeros (4,4), ( A2-Kf 2*GAM2 ' *Ce2 ) , . . 

Kf 2*GAM2 ' *Ce3 , zeros (4,6) ,Kf 2*GAM2 ’ *CeR] ; 

Atotal=[Atotal ;B3*Tl*Kcl ,B3*Tl*Kcl ,B3*T2*Kc2 ,B3*T2*Kc2 , . . 

( A3+B3STAR*Kc3 ) , B3STAR*Kc3 , zeros (6,8); . . 

Kf 3*GAM3'*Cel, zeros( 6, 6) ,Kf 3*GAM3 ' *Ce2 , zeros ( 6, 4) , zeros (6, 6) . . 

(A3-Kf3*GAM3'*Ce3) ,Kf 3*GAM3* *CeR] ; 

Atotal= [Atotal ;BR*Tl*Kcl,BR*Tl*Kcl,BR*T2*Kc2,BR*T2*Kc2, . . 
BR*T3*Kc3 ,BR*T3*Kc3 , AR] ; 

%SORT  Atotal  MATRIX  INTO  MODAL  SEQ  BY  ROWS  FIRST 
Ahold=zeros (48 , 40) ; 

POSCON=n(l:8,i) ; 

POSRES=n(9:12,i) ; 

RM=POSRES; 

MIXED= [ 1  2  3  13  14  21  22  23]; 

MIXEDR=  [ 33  34  35  36],* 

CTY  = [ 3  3  3  2  2  3  3  3]; 
for  h=l:l:4 

Ahold(POSCON , : )=Atotal (MIXED, : ); 

POSCON = POSCON  +  12; 

MIXED=MIXED+CTY; 

end 

for  g=l:l:2 

Aho 1 d ( POSRES , : )=Atotal (MIXEDR, : ) ; 

POSRES=POSRES+12; 

MIXEDR=MIXEDR+4 ; 
end 

%REMOVING  ROWS  CORRESPONDING  TO  NONEXISTANT  RESIDUAL  ERROR 
% STATES 

RM= [ ( RM+  2  4 ) , ( RM+  36)]; 

Ahold(RM, : )=[]; 

%SORT  Atotal  INTO  MODAL  SEQ  BY  COLUMNS  USING  Ahold 
Asort=zeros(40,48) ; 

POSCON=n(l:8,i) ; 

POSRES=n(9:12,i); 

RM= POSRES; 

MIXED= [ 1  2  3  13  14  21  22  23]; 

MIXEDR=[33  34  35  36]; 

CTY  = [ 3  3  3  2  2  3  3  3]; 
for  h=l:l:4 

Asort ( : , POSCON )=Ahold( : , MIXED) ; 

POSCON = POSCON  +  12; 

MIXED=MIXED+CTY; 

end 

for  g=l:l:2 

Asort ( : , POSRES )=Ahold( : , MIXEDR ) ; 

POSRES=POSRES+12; 

M I XEDR-M I XEDR + 4 ; 
end 

%REMOVING  COLUMNS  CORRESPONDING  TO  NONEXISTANT  RESIDUAL  ERROR 
%STATES 

RM=[ (RM+24) , (RM+36) ] ; 

Asort( : ,RM)=[ ] ; 

%REDUCING  MEMORY  REQUIREMENTS  BY  CLEARING  VARIABLES 
clear  W  ZW; 
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clear  ZR  ZC  K  L7  L8  L9  L10  Lll  L12; 

clear  K7  K8  K9  K10  Kll  K12  K7E  K8E  K9E  KlOE  K11E  K12E; 

clear  K7D  K8D  K9D  KIOD  K11D  K12D  Dn; 

clear  Atotal  C  Cel  Ce2  Ce3  CeR; 

clear  Ahold  POSCON  POSRES  RM  MIXED  CTY  MIXEDR; 

%SIMULATE  SYSTEM  RESPONSE  TO  AN  IMPULSE  AT  LEG10  (4th  ACTOR) 
%NOTE  DUE  TO  SYMMETRY  ONLY  AN  INPUT  AT  ANY  ONE  OF  THE  SIX 
%LEGS  PRODUCES  SIMILAR  OUTPUT  ONLY  IN  DIFFERENT  DIRECTIONS 
%  x  IS  THE  STATE  TIME  HISTORY,  y  IS  WHAT  THE  SENSORS  SEE 
%T  IS  THE  TIME  INCREMENT  AND  INTERVAL  SPECIFIED,  Dtotal  IS 
%THE  FEEDFORWARD  PORTION 
T=0 :  .02 : 20 ; 

Dtotal=[0  0000  0]*; 

[y,x]=impulse(Asort ,Btotal ,Ctotal , Dtotal ,4,T) ; 

%CONVERT  THE  TIME  HISTORY  STATE  MATRIX  X  TO  GENERALIZED  COORD 
%Q1  AND  Q2  ASSOCIATED  WITH  THE  LOS  OF  NODE  1 
COUNTER=l: length(T) ; 

WW=[1  00000000000]; 

WX=[0  1  0  0  0  0  0  0  0  0  0  0]; 

Ql=WW*U* ( X ( COUNTER ,(1:12)))'; 

Q2=WX*U* ( x( COUNTER ,(1:12)))'; 

Vl=WW*U*(x(COUNTER, (13:24) ) ) * ; 

V2=WX*U*(x(COUNTER, (13:24)) ) ' ; 

VELOCITY=sqrt (VI . “2  +  V2.*2); 

RADIUS=sqrt (Q1 . “2  +  Q2.*2); 

ANGLE=atan2(Q2,Ql) ; 

%INITI AL  SCREENING  CALCULATIONS,  FIRST  TWO  MOMENTS  OF 
%RAD  AND  ANGLE  ACROSS  TOTAL  TIME  (NOT  ERGODIC)  APPLIED  ONLY 
%WHEN  NO  NOISE  APPLIED 
%RADMEAN=mean( RADIUS) ; 

%RADSTD=std(RADIUS ) ; 

%FOLLOWING  CODE  HELPS  RELATE  SPATIAL  EVECTS  OF  VARIANCE 
%CON VERTING  POSITION  ANGLE  TO  BE  BETWEEN  0  AND  2PI 
%for  p=l : length (ANGLE) 

%  if  ANGLE(p)<0 
%  ANGLE(p)=ANGLE(p)  +  2*pi; 

%  end 
%end 

%ANGMEAN=mean ( ANGLE ) ; 

%ANGSTD=std (ANGLE) ; 

%% INITIALIZES  COUNTER  TO  ZERO 
%%  if  j==l 

%%  RADl=xeros(4 , length(RADIUS) ) ; 

%%  ANGl=zeros( 4 , length( ANGLE) ) ; 

%%  QUADl=zeros(4, length(ANGLE) ) ; 

%%  RADTOTAL=zeros(RADl ) ; 

%%  ANGTOTAL= zeros ( ANG1 ) ; 

%%  QUADTOTAL= zeros (QUAD1 ) ; 

%%  RADSQUARED=  z  er os ( RADI ) ; 

%  zero  out  the  tr(j)  used  for  each  noise  setting 

%%  end 

%%DIVIDING  THE  ANGLE  AND  RADIUS  INTO  QUADRANTS  1  THRU  4 
%%  for  k=l: LENGTH (ANGLE) 

%%  if  ANGLE(k)>3*pi/2 

%%  RADl(4,k)=RADIUS(k) ; 


%%  ANG1 ( 4 , k ) = ANGLE ( k ) ; 

%%  QUADl(4,k)=l; 

%%  else  if  ANGLE(k)>pi 
%%  RADl(3,k)=RADIUS(k) ; 

%%  ANG1 ( 3 , k)=ANGLE(k) ; 

%%  QUADl(3,k)=l; 

%%  else  if  ANGLE(k) >pi/2 

%%  RADI (2 , k)=RADIUS(k) ; 

%%  ANG1 ( 2 , k ) = ANGLE ( k ) ; 

%%  QUADl(2,k)=l; 

%%  else 

%%  RAD1(1 , k)=RADIUS(k) ; 

%%  ANG1 ( 1 , k ) = ANGLE ( k ) ; 

%%  QUAD1(1 ,k)=l ; 

%%  end 

%%  end 

%%  end 

%%  end 

%%R ADTOTAL = RADTOTAL + RADI ; 

%%RADSQUARED=RADSQUARED+RAD1  /2; 

%% AN  GTOT AL = AN  GTOT AL + AN G 1 ; 

%%QU  ADTOTAL  s QU ADTOTAL + QU  AD 1 ; 

%CALCULATE  A  MOVING  AVERAGE  FOR  EACH  CONRTOLLER  AT 
%EACH  NOISE  LEVEL.  STATISTIC  ARE  THEN  CALCULATED  FOR 
%THE  MOVING  AVERAGE 
f lag=0; 
if  flag==0 
for  k=927 : -X : 2 

ma(k)=mean(RADIUS( : , (k:k+74) ) ) ; 
if  flag==0  &  .05-ma(k)<= .002  ,%criteria 

% 

%IF  THE  NEW  Btotal  IS  USED  CHANGE  THE  ABOVE  .05  TO  SOME 
%NUMBER  BEWTEEN  .01  AND  .005  AS  LESS  ENERGY  IS  APPLIED 
% 

if  k<=926 

tr=(k+38)*.02;  ,%midpoint 
%  Va vg=mean ( VELOC I TY ( : ,(k:k+74))); 

else 
tr=25; 

%  Vavg=25; 

end 

f lag=l; 

end 

end 

end 

eval ( [ 'COMB' , int2str(i ) , ' (qc, j)=tr; ' ] ) 

%THIS  end  GOES  BACK  TO  THE  NOISE  LOOP  COUNTER  j 
end 

%THIS  end  GOES  BACK  TO  THE  Qc  LOOP  COUNTER  qc 
end 

%SAVES  SEPARATE  VARIABLE  FOR  EACH  COMBINATION 

eval ( [ ' save  ' , 'COMB1 , int2str(i) , '  COMB' , int2str( i ) ] ) 

% INITIALIZING  AN  OVERALL  VARIABLE  SCREEN210 
%TO  HOLD  ALL  OF  THE  INFORMATION  GENERATED 
%if  i==l 
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%  SCREEN210=zeros( 6 , 210 ) ; 
%end 


%eval ( [ ' SCREEN 2 10 ( 1 , 
%eval ( [ ' SCREEN210 ( 2 , 
%eval ( [ * SCREEN210 ( 3 , 
%eval ( [ ' SCREEN 2 10 ( 4 , 
%eval([ ’SCREEN210(5, 
%eval ( [ ' SCREEN 2 10 ( 6 , 


,  mt2str(  l ) ,  )  = 
,int2str(i ) , *  )  = 
,int2str(i) , ' )= 
, int2str(i ) , ' )= 
,int2str(i) , '  )  = 
,int2str(i) , ' )= 


Vavg  ;']) 
RADMEAN  ; * ] ) 
RADSTD 
ANGHEAN 
ANGSTD 


’]) 
;’]) 
’]) 

SO  AS 


%%SETTING  ANY  ZERO  ENTRIES  OF  QUADTOTAL  TO  1  SO  AS  TO 
%%AVOID  DIVIDING  BY  ZERO.  QUOTIENT  STILL  ZERO  SINCE  NUM  WILL 
%%BE  ZERO.  ALSO  MUST  AVOID  TAKING  THE  LOG  OF  ZERO 


%%for  p=l:4 

%%  for  k=l:length( QUADTOTAL) 

%%  if  QUADTOTAL ( p , k )  =  =  0 
%%  QUADTOTAL ( p , k ) = 1 ; 

%%  RAD_SN ( p , k ) =  0 ; 

%%  else 

%%  RAD_SN (p,k)=-10*log( RADSQU ARED ( p , k ) / QUADTOTAL ( p , k ) ) ; 

%%  end 
%%  end 
%%end 


%%FINDING  MEASURE  OF  PERFORMANCE 
%%ANG_MEAN=ANGTOTAL . /QUADTOTAL ; 

%%R AD_ME AN = RADTOT AL . /QUADTOTAL ; 

%%RAD_VAR=RADSQUARED . /QUADTOTAL  -RAD_MEAN . “ 2 ; 

%%eval ( [ ’ save  ' , 'RAD_MEAN' ,int2str(i) , '  RAD_MEAN ' ] ) ; 

%%eval([ 'save  ' , 'RAD_VAR' , int2str( i ) , '  RAD_VAR ' ] ) ; 

%%eval ( [ ' save  ' , 'RAD_SN' , int2str(i) , '  RAD_SN']); 

%%eval ( [ ' save  ' , ' ANG_MEAN' , int2str(i ) , '  ANG_MEAN']); 

%clear  RAD_MEAN  RAD_VAR  RAD_SN  ANG_MEAN 

%clear  ANGTOTAL  QUADTOTAL  RADTOTAL  RADSQUARED  RADI  ANG1  QUAD1 

%clear  Ctotal  Asort  x  y  Q1  Q2  VI  V2  ANGLE  RADIUS 

%%THIS  end  GOES  BACK  TO  SELECT  ANOTHER  CONTROLLER  ARRANGEMENT 

end 

%save  SCREEN210  SCREEN 2 10 
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Appendix  G 
Regression  Results 


0-1 


Table  G-l 


SCREENING  DESIGN  FOR  Q  ,  and  Q _  ARR  126-47-358 

or  oc 


GROUP  1 

STATE 

WEIGHTING 

COEFF . 

CONTROLLER  GAIN 

GROUP  2  £2 

STATE 

WEIGHTING 

COEFF . 

ESTIMATOR  GAIN 

• 

GROUP  l2 

hi 

GROUP  12 

^12 

GROUP  2 2 

l22 

• 

GROUP  122 

hl2 

GROUP  12 2 

5122 

GROUP  1222 

?1122 

USED  32  FULL  FACTORIAL 
SIGNIFICANT  GROUPS  1,  2  AND  12 
LACK  OF  FIT 


9  RUNS  X  NOISE 
(F=159 , 120 , 37 ) 
(F=1.4(10_26)) 


EXPERIMENTAL  DESIGN  MATRIX  Z 


(3522-2)x 

(3f11-2)52  (3e22-2)?1  (3^11-2) 


-1  -1  1 

2  0-2 

-111 
0  2-2 

0  0  4 

0  -2  -2 

1-11 
-2  0  -2 

111 
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Table  G-2 


ANOVA  FOR  Q  -  AND  <5  SCREENING 
01  oc 

MEAN  RESPONSE  ARR  126-47-358 


MEAN 

GROUP  1 

GROUP  2 

GROUP  l2 

ss 

6 . 7030e+03 

2 . 1942e+01 

1 . 6653e+01 

1 . 1574e-01 

df 

1 . 0000e+00 

1 . 0000e+00 

1 . 0000e+00 

1.0000e+00 

MS 

6 . 7030e+03 

2 . 1942e+01 

1 . 6653e+01 

1 . 1574e-01 

F 

4 . 8628e+04 

1 . 5918e+02 

1 . 2081e+02 

8 . 3965e-01 

0 

5 • 2521e+00 

-3 . 6802e-01 

-3 . 2062e-01 

1 . 5432e-02 

GROUP  12 

GROUP  2 2 

GROUP  122 

GROUP  122 

SS 

5 . 0700e+00 

8 . 4807e-02 

1 . 2346e+00 

1 . 9753e-03 

df 

1 . 0000e+00 

1 . 0000e+00 

1 . 0000e+00 

1 . 0000e+00 

MS 

5 . 0700e+00 

8 . 4807e-02 

1 . 2346e+00 

1.9753e-03 

F 

3 . 6781e+01 

6 . 1524e-01 

8 . 9563e+00 

1 . 4330e-02 

0 

-2 . 1667e-01 

-1 . 3210e-02 

6 . 1728e-02 

2 . 4691e-03 

GROUP  1222 

LACK  FIT 

PURE  ERROR 

TOTAL 

SS 

2 . 8624e-01 

1. 9382e-27 

3 . 2255e+01 

6 . 7807e+03 

df 

1 . 0000e+00 

1.0000e+00 

2 . 3400e+02 

2 . 4300e+02 

MS 

2 . 8624e-01 

1 . 9382e-27 

1 . 3784e-01 

2 . 7904e+01 

F 

2 . 0765e+00 

1 . 4061e-26 

0 

0 

0 

1 . 7160e-02 

0 

0 

5 . 2521e+00 

SIGNAL-TO-NOISE  ARR  126-47-358 


MEAN 

GROUP  1 

GROUP  2 

GROUP  l2 

• 

SS 

4 . 1609e+03 

5 . 7928e-01 

4 . 5764e-01 

2 . 0349e-03 

df 

1 . 0000e+00 

1 . 0000e+00 

1 . 0000e+00 

1 . 0000e+00 

MS 

4 . 1609e+03 

5.7928e-01 

4 . 5764e-01 

2 . 0349e-03 

F 

not  applicable:  9  coefficients  for  9 

data  points 

0 

-2 . 1502e+01 

3 . 1072e-01 

2.7618e-01 

-1.0633e-02 

• 

GROUP  12 

GROUP  2 2 

GROUP  122 

GROUP  122 

SS 

1 . 5193e-01 

4 . 0050e-03 

2 . 8526e-02 

1 . 5059e-05 

df 

1 . 0000e+00 

1 . 0000e+00 

1.0000e+00 

1 . 0000e+00 

MS 

1 . 5193e-01 

4.0050e-03 

2 . 8526e-02 

1 . 5059e-05 

F 

** 

** 

** 

** 

0 

1 . 9489e-01 

1 . 4916e-02 

-4 . 8756e-02 

1 . 1202e-03 

• 

GROUP  1222 

LACK  FIT 

TOTAL 

SS 

7 . 9895e-03 

** 

4 . 1621e+03 

df 

1.0000e+00 

** 

9.0000e+00 

MS 

7 . 9895e-03 

** 

4 . 6245e+02 

F 

** 

** 

0 

• 

0 

-1 . 4897e-02 

** 

-2 . 1502e+01 

G-3 


Table  G-3 


SCREENING  DESIGN  POR  and  R  ARR  126-47-358 

oc  oc 

GROUP  1  STATE  WEIGHTING  COEFF.  CONTROLLER  GAIN 

GROUP  2  £2  CONTROL  WEIGHTING  COEFF.  CONTROLLER  GAIN 

GROUP  l2  5X1 

GROUP  12  512 

GROUP  22  l22 

GROUP  122  ?112 

GROUP  122  ?122 

GROUP  1222  €1122 

USED  3 2  FULL -ORDER  9  RUNS  X  NOISE 

SIGNIFICANT  GROUPS  1,  2  AND  12  (F=1371 , 1407 , 345) 

LACK  OF  FIT  (F=2 . 3(10_25) ) 

EXPERIMENTAL  DESIGN  MATRIX  Z 

(3£22-2)x 


«1 

e2 

3?ir2 

512 

H22~2 

(3^11-2)52 

(35n-2) 

-1 

-1 

1 

1 

1 

-1 

-1 

1 

0 

-1 

-2 

0 

1 

2 

0 

-2 

1 

-1 

1 

-1 

1 

-1 

1 

1 

-1 

0 

1 

0 

-2 

0 

2 

-2 

0 

0 

-2 

0 

-2 

0 

0 

4 

1 

0 

1 

0 

-2 

0 

-2 

-2 

-1 

1 

1 

-1 

1 

1 

-1 

1 

0 

1 

-2 

0 

1 

-2 

0 

-2 

1 

1 

1 

1 

1 

1 

1 

1 

G-4 


Table  G-4 


ANOVA  FOR  O  „  AND  R  SCREENING 
oc  oc 

MEAN  RESPONSE  ARR  126-47-358 


MEAN 

GROUP  1 

GROUP  2 

GROUP  l2 

ss 

3 . 3169e+04 

2 . 0299e+02 

2 . 0826e+02 

3 . 0658e-02 

df 

1.0000e+00 

1 . 0000e+00 

1 . 0000e+00 

1 . 0000e+00 

MS 

3 . 3169e+04 

2 . 0299e+02 

2 . 0826e+02 

3 . 0658e-02 

F 

2 . 2407e+05 

1 . 3713e+03 

1 . 4069e+03 

2.0710e-01 

0 

1 . 1683e+01 

-1 . 1194e+00 

1 . 1338e+00 

7 . 9424e-03 

GROUP  12 

GROUP  2 2 

GROUP  122 

GROUP  122 

SS 

5 . 1143e+01 

3 . 2253e-01 

1 . 9497e+01 

1 . 7269e+01 

df 

1 . 0000e+00 

1 .0000e+00 

1 . 0000e+00 

1 . 0000e+00 

MS 

5 . 1143e+01 

3 . 2253e-01 

1 . 9497e+01 

1 . 7269e+01 

F 

3 . 4549e+02 

2 . 1788e+00 

1 . 3171e+02 

1 . 1666e+02 

0 

6.8815e-01 

-2 . 5761e-02 

-2 . 4531e-01 

2 . 3086e-01 

GROUP  1222 

LACK  FIT 

PURE  ERROR 

TOTAL 

SS 

4 . 3506e+01 

3 . 3312e-26 

3 . 4639e+01 

3 . 3747e+04 

df 

1.0000e+00 

1 . 0000e+00 

2 . 3400e+02 

2 . 4300e+02 

MS 

4 . 3506e+01 

3 . 3312e-26 

1 . 4803e-01 

1 . 3888e+02 

F 

2 . 9390e+02 

2 . 2504e-25 

0 

0 

0 

2  •  1156e-01 

0 

0 

1 . 1683e+01 

SIGNAL-TO-NOISE  ARR  126-47-358 


MEAN 

GROUP  1 

GROUP  2 

GROUP  l2 

• 

SS 

5 . 6030e+03 

1 . 1375e+00 

1 . 1774e+00 

1 . 0033e-04 

df 

1 . 0000e+00 

1 . 0000e+00 

1.0000e+00 

l.OOCOe+OO 

MS 

5 . 6030e+03 

1 . 1375e+00 

1 . 1774e+00 

1 . 0033e-04 

F 

not  applicable:  9  coefficients  for  9 

data  points 

0 

-2 . 4951e+01 

4 . 3542e-01 

-4 . 4298e-01 

-2 . 3609e-03 

• 

GROUP  12 

GROUP  2 2 

GROUP  122 

GROUP  122 

SS 

3 . 1281e-01 

2 . 4659e-03 

1 . 0867e-01 

9. 2076e-02 

df 

1 . 0000e+00 

1 . 0000e+00 

1.0000e+00 

1 . 0000e+00 

MS 

3 . 1281e-01 

2. 4659e-03 

1 . 0867e-01 

9. 2076e-02 

F 

** 

** 

** 

** 

0 

-2 . 7965e-01 

1 . 1705e-02 

9 . 5163e-02 

-8 . 7596e-02 

• 

GROUP  1222 

LACK  FIT 

TOTAL 

SS 

2 . 5324e-01 

** 

5 . 6061e+03 

df 

1.0000e+00 

** 

9.0000e+00 

MS 

2 . 5324e-01 

*  * 

6 . 2290e+02 

F 

** 

** 

0 

• 

0 

-8 . 3872e-02 

** 

-2 . 4951e+01 

0-5 


Table  0-5 


•  GROUP  SCREEN  DESIGN  FOR  CONTROLLER  GAIN 


GROUP 

1 

*1 

SUBCONTROLLER 

1 

POSITION 

3 

MODES 

GROUP 

2 

*2 

SUBCONTROLLER 

1 

VELOCITY 

3 

MODES 

• 

GROUP 

3 

*3 

SUBCONTROLLER 

2 

POSITION 

2 

MODES 

GROUP 

4 

*4 

SUBCONTROLLER 

2 

VELOCITY 

2 

MODES 

GROUP 

5 

^5 

SUBCONTROLLER 

3 

POSITION 

3 

MODES 

• 

GROUP 

6 

*6 

SUBCONTROLLER 

3 

VELOCITY 

3 

MODES 

USED  2?“2  FRACTIONAL  3x16  RUNS  X  NOISE 

BEST  COMBINATION  126-45-358 

SIGNIFICANT  GROUPS  2  AND  6  (F=10,110) 

LACK  OF  FIT  (F=.88) 

EXPERIMENTAL  DESIGN  MATRIX  Z 


*5  *6 

1  -1 

1  -1 

1  1 

1  1 

1  1 

1  1 

1  -1 

1  -1 

1  1 

1  1 

1  -1 

1  -1 

1  *  j. 

1  -1 

1  1 

1  1 


G-6 


Table  G-6 


ANOVA  FOR  CONTROLLER  GAIN  GROUP  SCREEN 
MEAN  RESPONSE  ARR  126-47-358 


MEAN 

GROUP  1 

GROUP  2 

GROUP  3 

ss 

df 

MS 

F 

0 

1 . 3343e+04 

1 . 0000e+00 

1 . 3343e+04 

2 . 7478e+05 

5 . 5576e+00 

2 • 9668e-02 

1 . 0000e  +  00 

2 . 9668e-02 

6 . 1094e-01 

8 . 2870e-03 

4 . 8669e-01 
1.0000e+00 

4 . 8669e-01 

1 . 0022e+01 
-3 . 3565e-02 

2 . 0833e-04 
1 . 0000e+00 
2.0833e-04 
4 . 2901e-03 
-6 . 9444e-04 

GROUP  4 

GROUP  5 

GROUP  6 

ss 

df 

MS 

F 

0 

4 . 1565e-03 

1 . 0000e+00 

4 . 1565e-03 
8.5593e-02 
-3 . 1019e-03 

1 . 2245e-02 

1 . 0000e+00 

1 . 2245e-02 

2 . 5217e-01 
-5 . 3241e-03 

5 . 3734e+00 
1.0000e+00 
5 . 3734e+00 
1 . 1065e+02 
-1 . 1153e-01 

LACK  FIT 

PURE  ERROR 

TOTAL 

SS 

df 

MS 

F 

0 

3.8628e-01 

9 . 0000e+00 
4.2920e-02 

8 . 8384e-01 

0 

2 . 0201e+01 

4 . 1600e+02 

4 . 8561e-02 

0 

0 

1 . 3370e+04 
4 . 3200e+02 
3 . 0949e+01 
0 

5 . 557  6e+00 

SIGNAL-TO-NOISE 

ARR  126-47-358 

MEAN 

GROUP  1 

GROUP  2 

GROUP  3 

SS 

df 

MS 

F 

0 

1 . 3657e+04 

1 . 0000e+00 

1 . 3657e+04 

3 . 5846e+06 
-2 . 9216e+01 

2 ,0579e-03 
1 . 0000e+00 
2 . 0579e-03 
5 . 4013e~01 
-1 . 1341e-02 

4 . 4054e-02 
1.0000e+00 

4 . 4054e-02 

1 . 1563e+01 

5 . 2473e-02 

4 . 0439e-05 
1 . 0000e+00 
4 . 0439e-05 
1 . 0614e-02 
1 . 5898e-03 

GfvOUP  4 

GROUP  5 

GROUP  6 

SS 

df 

MS 

F 

0 

3 . 0570e-04 
1.0000e+00 
3 . 0570e-04 
8 . 0236e-02 
4 . 3711e-03 

1 . 2396e-03 

1 . 0000e+00 

1 . 2396e-03 

3 . 2534e-01 

8 . 8018e-03 

4 . 5181e-01 
1 . 0000e+00 
4 . 5181e-01 
1 . 1859e+02 
1 . 6804e-01 

LACK  FIT 

TOTAL 

SS 

df 

MS 

F 

0 

3 . 4290e-02 
9.0000e+00 

3 . 8100e-03 

0 

0 

1 . 3658e+04 
1 . 6000e+01 
8.5362e+02 
0 

-2 . 9216e+01 

G-7 


Table  G-7 


INDIVIDUAL  SCREEN  DESIGN  FOR  CONTROLLER  GAIN 
(  1  ARRANGEMENT,  SET  GROUPS  TO  MIN  VARIANCE  ) 

PACTOR  1  SUBCONTROLLER  1  VELOCITY  MODE  1 

PACTOR  2  f2  SUBCONTROLLER  1  VELOCITY  MODE  2 

PACTOR  3  £3  SUBCONTROLLER  1  VELOCITY  MODE  6 

PACTOR  4  £4  SUBCONTROLLER  3  VELOCITY  MODE  3 

PACTOR  5  £5  SUBCONTROLLER  3  VELOCITY  MODE  5 

PACTOR  6  SUBCONTROLLER  3  VELOCITY  MODE  8 

USED  2?~2  FRACTIONAL  16  RUNS  X  NOISE 

WITH  CENTER  POINTS 

SIGNIFICANT  FACTORS  4  AND  6  (F=50,56) 

NOT  SIGNIFICANT,  BUT  INCLUDED  FACTOR  2  (F=1.6) 

LACK  OF  FIT  CROSS  TERMS  (F  1Q  -  1.63)  (F=3.4) 

LACK  OF  FIT  QUADRATIC  TERMS  (F  1Q  *  2.71)  (F=1.8) 

EXPERIMENTAL  DESIGN  MATRIX  Z 

«1  *2  *3  *4  *5  e6 

-1  -1  -1  -1  -1  -1 

1  -1  -1  -1  1  -1 

-1  1-1-1  1  1 

1  1  -1  -1  -1  1 

-1-1  1-1  1  1 

1-1  1-1-1  1 
-1  1  1-1  -1  -1 
1  1  1-1  1-1 
-1  -1-1  1-1  1 

1-1-1  1  1  1 

-1  1-1  1  1-1 

1  1-1  1-1-1 

-1-1  1  1  1-1 

1-1  1  1-1-1 

-1111-11 
111111 
0  0  0  0  0  0 


G-8 


Table  G-8 


ANOVA  FOR  CONTROLLER  GAIN  INDIVIDUAL  SCREEN 
MEAN  RESPONSE  ARR  126-47-358  . 

CONTROLLER  1  VELOCITY 


MEAN 

MODE  1 

MODE  2 

MODE  6 

ss 

1 . 3392e+04 

1.0800e-02 

5 . 1570e-02 

8 . 1815e-03 

df 

1.0000e+00 

1 . 0000e+00 

1 . 0000e+00 

1 . 0000e+00 

MS 

1 . 3392e+04 

1 . 0800e-02 

5 . 1570e-02 

8 . 1815e-03 

F 

4 . 1637e+05 

3 . 3578e-01 

1 . 6033e+00 

2 . 5437e-01 

0 

5 . 5678e+00 

5 . 0000e-03 

-1.0926e-02 

-4 . 3519e-03 

CONTROLLER 

3  VELOCITY 

MODE  3 

MODE  5 

MODE  8 

SS 

1 . 6084e+00 

3 .0000e-04 

1 . 8148e+00 

df 

1 . 0000e+00 

1 . 0000e+00 

1 . 0000e+00 

MS 

1 . 6084e+00 

3 . 0000e-04 

1 . 8148e+00 

F 

5 . 0007e+01 

9 . 3271e-03 

5 . 6423e+01 

0 

-6 . 1019e-02 

-8 . 3333e-04 

-6 . 4815e-02 

LACK 

OF  FIT 

CROSS  PROD 

PURE  QUAD 

PURE  ERROR 

TOTAL 

SS 

9 • 9502e-01 

5 . 8008e-02 

1 . 4217e+01 

1 . 3411e+04 

df 

9 . 0000e+00 

1.0000e+00 

4 . 4200e+02 

4 . 5900e+02 

MS 

1 . 1056e-01 

5 . 8008e-02 

3 . 2164e-02 

2 . 9218e+01 

F 

3 . 4373e+00 

1 . 8035e+00 

0 

0 

0 

0 

0 

0 

5 . 5678e+00 

SIGNAL-TO-NOISE 

1  ARR  126-47-358 

CONTROLLER 

1  VELOCITY 

MEAN 

MODE  1 

MODE  2 

MODE  6 

SS 

1 . 3671e+04 

7. 9481e-04 

5 . 1101e-03 

9. 8009e-04 

df 

1 . 0000e+00 

1 . 0000e+00 

1 . 0000e+00 

1.0000e+00 

MS 

1 . 3671e+04 

7 . 9481e-04 

5 . 1101e-03 

9 . 8009e-04 

F 

1 . 4216e+06 

8 . 2649e-02 

5 . 3139e-01 

1.0192e-01 

0 

-2 . 9231e+01 

-7 . 0481e-03 

1.7871e-02 

7 . 8266e-03 

CONTROLLER 

3  VELOCITY 

MODE  3 

MODE  5 

MODE  8 

SS 

1 . 3469e-01 

8 . 9692e-05 

1 . 5108e-01 

df 

1.0000e+00 

1.0000e+00 

1 . 0000e+00 

MS 

1 . 3469e-01 

8 . 9692e-05 

1 . 5108e-01 

F 

1 . 4006e+01 

9 . 3267e-03 

1 . 5710e+01 

0 

9 . 1752e-02 

2 . 367  6e-03 

9.7173e-02 

LACK  FIT 

TOTAL 

SS 

8 . 6550e-02 

1 . 3671e+04 

df 

9.0000e+00 

1 . 6000e+01 

MS 

9 . 6166e-03 

8 . 5445e+02 

F 

0 

0 

0 

0 

-2 . 9231e+01 

G-9 


Table  G-9 


CONTROLLER  GAIN  INTERACTION  DESIGN 
(OTHER  PACTOR  SET  TO  MINIMIZE  VARIANCE) 


FACTOR  1 
PACTOR  2 
PACTOR  3 
FACTOR  12 
FACTOR  13 
FACTOR  23 


*1 

*2 

*3 

512 

513 

*23 


SUBCONTROLLER  1  VELOCITY  MODE  2 
SUBCONTROLLER  3  VELOCITY  MODE  3 
SUBCONTROLLER  3  VELOCITY  MODE  8 


8  RUNS  x  NOISE 
(F=29, 20,14) 
(F= . 64) 


USED  PULL  FACTORIAL 
SIGNIFICANT  TERMS  2  3  AND  23 
LACK  OF  FIT 

TIMER  =  5.5669  -.05676Qy3  -.0475Qv8  -.0397Qv3Qv8 
EXPERIMENTAL  DESIGN  MATRIX  Z 
^1  £3  ^13  ^23 


-1 

1 

-1 

1 

-1 

1 

-1 

1 


-1 

-1 

1 

1 

-1 

-1 

1 

1 


-1 

-1 

-1 

-1 

1 

1 

1 

1 


1 

-1 

-1 

1 

1 

-1 

-1 

1 


1 

-1 

1 

-1 

-1 

1 

-1 

1 


1 

1 

-1 

-1 

-1 

-1 

1 

1 


G-10 


Table  G-10 


ANOVA  POR  CONTROLLER  GAIN  INTERACTION 
MEAN  RESPONSE  ARR  126-47-358 


MEAN 

MODE  ? 

MODE  3 

MODE  8 

ss 

6 . 6940e+03 

4 . 0017e-02 

6 . 9587e-01 

4.8'»35e-01 

df 

1.0000e+00 

1.0000e+00 

1 . 0000e+00 

1.0000e+00 

MS 

6 . 6940e+03 

4.0017e-02 

6 . 9587e-01 

4 . 8735e-01 

F 

2 . 7969e+05 

1 . 6720e+00 

2 . 9075e+01 

.  2.0363e+01 

e 

5 . 5669e+00 

-1 . 3611e-02 

-5 . 6759e-02 

-4 . 7500e-02 

CROSS  2-3 

CROSS  2-8 

CROSS  3-8 

SS 

1 . 6713e-02 

1 . 8150e-02 

3 . 4082e-01 

df 

1 . 0000e+00 

1 . 0000e+00 

1.0000e+00 

MS 

1 . 6713e-02 

1 . 8150e-02 

3 . 4082e-01 

P 

6 . 9831e-01 

7 . 5836e-01 

1 . 4240e+01 

e 

-8 . 7963e-03 

-9 . 1667e-03 

-3 . 9722e-02 

LACK  FIT 

PURE  ERROR 

TOTAL 

SS 

1 . 5335e-02 

4. 9781e+00 

6 . 7006e+03 

df 

1 .0000e+00 

2 . 0800e+02 

2 . 1600e+02 

MS 

1 . 5335e-02 

2 . 3933e-02 

3 . 1021e+01 

F 

0 

0 

0 

e 

0 

0 

5 . 5669e+00 

SIGNAL-TO  NOISE 

ARR  126-47- 

358 

MEAN 

MODE  2 

MODE  3 

MODE  8 

ss 

6 . 8344e+03 

3. 6282e-03 

5 . 8972e-02 

4 . 0888e-02 

df 

1.0000e+00 

1.0000e+00 

1 . 0000e+00 

1.0000e+00 

MS 

6 . 8344e+03 

3 . 6282e-03 

5 . 8972e-02 

4 . 0888e-02 

p 

4 . 8434e+06 

2 . 5712e+00 

4 . 1793e+01 

2 . 8977e+01 

e 

-2 . 9228e+01 

2 . 1296e-02 

8 . 5858e-02 

7 . 1492e-02 

CROSS  2-3 

CROSS  2-8 

CROSS  3-8 

SS 

1 . 5391e-03 

1 . 6650e-03 

2 . 8344e-02 

df 

1.0000e+00 

1.0000e+00 

1.0000e+00 

MS 

1 . 5391e-03 

1 . 6650e-03 

2 . 8344e-02 

F 

1 . 0908e+00 

1 . 1799e+00 

2.0087e+01 

e 

1 . 3871e-02 

1 . 4426e-02 

5 . 9523e-02 

LACK  PIT 

TOTAL 

ss 

1 . 4111e-03 

6 . 8345e+03 

df 

1.0000e+00 

8 . 0000e+00 

MS 

1.4111e-03 

8 . 5431e+02 

p 

0 

0 

0 

0 

-2 . 9228e+01 

0-11 
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Block  19 .( continued)  Abstract 


This  investigation  developed  and  tested  statistical  design 
methods  for  configuring  a  robust  decoupled  controller  on  a 
lightly  damped  structure.  The  test  model  used  was  a  lumped  mass 
finite  element  representation  of  the  Charles  Stark  Draper 
Laboratory  Model  I  (CSDL-I).  The  decoupled  control  methods 
consisted  of  a  system  of  three  individual  subcontrollers  designed 
on  the  basis  of  a  subset  of  plant  dynamics  with  control  authority 
enacted  through  the  total  system  of  structural  actuators. 
Transforms  were  applied  to  the  subcontrollers  to  insure  that 
dynamic  coupling,  called  observation  spillover  and  control 
spillover,  didn't  destabilize  the  global  system.  The  decoupled 
controller  system  was  designed  based  on  only  the  first  eight  of 
the  twelve  natural  frequencies  (modes)  for  the  CSDL-I.  The 
remaining  four  high  frequency  modes  were  modeled  as  residuals. 

This  investigation  used  signal-to-noise  ratios,  orthogonal 
numerical  experiments,  and  linear  regression  to  efficiently  probe 
the  design  space  and  to  produce  a  robust  control  system.  The 
measure  of  system  performance  was  the  structural  alignment  of  a 
particular  node  referred  to  as  1 ine-of -sight  (LOS).  This  LOS 
requirement  simulated  space  vehicle  requirements  associated  with 
aligning  antennas,  LASER  beams,  solar  arrays,  etc. 

It  was  concluded  that  the  techniques  were  efficient  in 
configuring  a  decoupled  controller  to  be  robust  to  sensor  noise 
and  modal  truncation  while  still  maintaining  adequate  LOS 
performance.  The  techniques  reduced  the  CPU  time  required  from  a 
possible  76,354  years  for  a  global  search  to  only  16.5  hours.  It 
was  also  found  that  the  robust  controller  used  a  significant 
amount  of  actuator  authority  and  it  was  recommended  that  further 
work  be  done  to  expand  the  statistical  design  techniques  to 
include  actuator  authority  as  an  output  metric. 
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